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This  document  reports  results  obtained  by  the  Lockheed-Georgia  Co.  and 
Grumman  Aerospace  Corp.  under  AFFDL  Contract  #  F33615-78-C-3014.  The  purpose  of 
the  contract  was  to  develop  and  validate  a  new  transonic  wing  design  procedure 
using  the  numerical  optimization  technique.  The  new  procedure  was  used  to  design 
both  a  transport  and  a  fighter  configuration. 

Because  the  missions  and  design  requirements  of  a  fighter  and  transport  are 
so  different,  the  design  procedure  was  developed  along  parallel  lines.  Lockheed- 
Georgia  Co.  developed  the  transport  design  procedure,  and  Grumman  Aerospace  Corp. 
did  the  fighter  design. 

This  is  Part  2  of  a  two-part  volume:  Part  1  details  the  transport  design  code 
Users  Guide  and  Part  2  is  the  fighter  design  code  Users  Guide.  There  are  two  other 
volumes  which  make  up  the  final  report.  Volume  1  is  an  executive  summary.  Volume 
1  is  also  divided  into  two  parts  with  Part  1  dealing  with  the  transport  design  and 
Part  2  concerned  with  the  fighter  design.  Volume  2  is  a  detailed  discussion  of 
the  technical  tasks  performed  under  this  contract.  Volume  2  is  similarly  divided 
into  two  parts. 

Personnel  who  contributed  to  this  contract  effort  are:  Lockheed-Georgia 
Company,  A.  J.  Srokowski,  M.  E.  Lores,  R.  A.  Weed  and  P.  R.  Smith;  Grumman 
Aerospace  Corp.,  P.  Aidala.  The  authors  also  wish  to  acknowledge  the  assistance 
given  by  Capt.  R.  A.  Large  who  was  the  AFWAL  contract  monitor. 
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SECTION  I 


INTRODUCTION 

Computational  aerodynamic  methods  for  the  design  of  future  configurations 
are  becoming  more  important.  New  technology  concepts  and  new  configurational 
concepts  are  being  developed  to  improve  transonic  performance.  Accurate  per¬ 
formance  predictions  are  needed  to  confidently  evaluate  the  benefits  of  new 
concepts.  In  addition,  the  pursuit  of  increased  performance  and  reduced  develop¬ 
ment  time  and  costs  intensify  the  use  of  three-dimensional  computational  analysis. 

The  fighter  design  code  developed  in  this  study  is  a  wing-body-canard 
transonic  analysis  capability  coupled  with  numerical  optimization.  The  multiple 
lifting  surface  transonic  analysis  capability  is  the  first  available  for  general 
applications.  The  wing-body- canard  analysis  capability  was  developed  from  the 
wing-body  code  of  Boppe  [l],  while  the  optimization  routines  are  essentially 
direct  copies  of  those  of  Vanderplaats  [8,9]. 

The  code  has  been  named  PANDORA — Preliminary  Automated  Numerical  Design  Of 
Realistic  Aircraft.  It  was  developed  on  a  CDC  7600  computer  at  NASA  Ames  Research 
Center.  Computer  time  for  the  code  development  and  application  was  provided 
through  the  Applied  Computational  Aerodynamics  (ACA)  Branch  at  Ames.  In  addition, 
parts  of  the  work  were  performed  at  Ames  through  the  Grumman-Ames  Research 
Associate  Program.  The  code  is  stored  on  the  Ames  central  facility  7600  in  CDC 
UPDATE  format.  Technical  details  of  accessing  the  code  can  be  provided  by  the  ACA 
Branch  at  Ames. 

This  Users  Guide  will  provide  the  information  needed  to  use  the  PANDORA  code 
for  configuration  analysis  or  aerodynamic  optimization.  Volume  2,  Part  2  of  this 
final  report  contains  important  discussions  of  the  code  accuracy  and  application. 
Unnecessary  details  of  the  code  which  are  the  same  as  described  in  References  [l], 
[8]  or  [9]  are  not  repeated  here.  Those  references  should  be  consulted  if  the 
particular  code  details  are  of  interest. 

SECTION  II 

WING- BODY- CANARD  TRANSONIC  ANALYSIS 

The  basis  computational  method  employed  in  the  transonic  analysis  is  that  of 
Boppe  [ 1 ] .  An  earlier  version  (1978)  of  the  wing-body  analysis  code  was  used  as 
the  base  on  which  to  develop  the  wing-body-canard  analysis  capability.  The  new 
developments  are  described  in  detail  here.  The  common  parts  (flow  equation, 
embedded  grid  interfacing,  body  modeling,  viscous  effects)  are  described  briefly. 
Reference  [1]  should  be  consulted  if  more  details  are  of  interest. 

1.  FLOW  EQUATION 

The  flow  equation  used  in  the  analysis  is  an  ’’extended"  small-disturbance 
equation: 

(■->£  -  <V1»6X  -  *T  »5>Xx  -  “fryV 

2  (,) 
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The  additional  terms  have  been  added  to  better  capture  swept  shock  waves  and  more 
accurately  determine  the  critical  velocity.  Empirical  modifications  and 
similarity  variables  are  not  employed. 

Pressure  coefficients  on  wing  surfaces  are  computed  using  the  following 
equation: 

Cp  =  '  (2*x  +  (1  '  +  ^  (2) 

To  simplify  velocity  computations  on  the  non-planar  body  surface,  a  simplified 
equation  is  used: 

cp  =  -  (2<t>x  +  (l-M2)*2)  (3) 

The  computational  space  used  in  the  present  method  is  illustrated  in  Fig¬ 
ure  1.  This  space  is  filled  with  a  relatively  crude  Cartesian  mesh.  Instead  of 
adopting  a  far-field  solution  for  the  grid  outer  boundaries,  the  original  x,  y,  z 
region  is  stretched  to  £,r),£  region  in  which  the  boundaries  correspond  to 
infinity.  The  flow  field  potential  is  set  to  zero  on  all  bounding  planes  except 
the  downstream  plane  for  which  the  following  equation  is  solved. 

<(>  +  4>  =  0  (4) 

Yyy  rzz 

The  following  conditions  are  enforced  at  the  synmetry  plane. 

♦  =  0  (5a) 

<()  =0  (5b) 

Yxy 

The  numerical  solution  of  the  flow  equation  uses  successive  line  over¬ 
relaxation  (SLOR) .  Vertical  line  solutions  proceed  streamwise  at  a  given  spanwise 
plane.  For  the  crude  grid  solution,  the  spanwise  planes  are  solved  sequentially 
starting  at  the  most  outboard  station.  For  the  embedded  fine  grid  solution,  the 
spanwise  grid  systems  are  solved  sequentially  from  the  wing  root  to  the  wing  tip 
and  then  from  the  canard  root  to  the  canard  tip. 

2.  GRID  GENERATION 

The  key  item  to  enable  wing-canard  analysis  without  undue  computer  storage 
requirements  was  the  development  of  suitable  grid  point  distributions.  The  global 
crude  mesh  uses  51,  26  and  31  planes  in  the  streamwise,  spanwise  and  vertical 
directions,  respectively.  This  results  in  41,106  mesh  points. 

In  order  to  distribute  the  crude  grid  streamwise  mesh  planes  as  effectively 
as  possible,  separate  grids  are  used  for  wing  and  wing-canard  analysis.  Without  a 
canard,  the  streamwise  grid  is  generated  with  a  tangent  function: 

X  =  A1TAN((n/2)4)  (6) 

This  transforms  the  finite  equally-spaced  computational  domain  -1  <  £  <  1  to  the 
physical  region  -»  <  X  <  ».  The  grid  is  centered  on  the  wing  planform. 
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When  a  canard  is  present,  the  grid  of  equation  (6)  usually  results  in  too  few 
crude  grid  planes  intercepting  the  canard  planform.  The  streamwise  grid 
transformation  used  for  wing-canard  combinations  is: 

x  =  a0  +  ai4  +  a/  +  a3C3  +  +  a5C5  +  as£/U-£2)  (7a) 

i  =  i/i^  (7b) 

This  transforms  the  finite  domain  -1  <  £  <  1  to  the  infinite  region  -»  <  X  < 
00 .  The  constant  a  is  empirical,  controlling  the  rate  of  stretching  near 

infinity.  The  points  at  ^  =  i  1  are  made  the  points  of  maximum  density  of  mesh 

planes  by  specifying  the  second  derivative  of  X  to  be  zero.  The  value  of  X  and  the 
first  derivative  of  X  at  £  =  ±  1  are  also  specified  to  determine  the  coefficients 
a^  to  a,..  The  value  of  is  adjusted  iteratively  to  place  approximately  fifty 
percent  of  the  total  mesh  planes  between  the  most  forward  and  the  most  aft  points 
on  the  wing-canard  combination. 

For  an  aft  swept  wing-canard  combination,  the  wing  tip  and  canard  tip 
determine  the  transformation.  The  mid-points  of  the  two  tip  chords  are  the  values 
of  X  for  £  =  ±  1.  The  first  derivative  of  X  is  set  to  result  in  a  nominal  six  mesh 
planes  intercepting  each  of  the  tip  chords.  An  example  of  the  resulting  wing- 
canard  grid  transformation  and  the  corresponding  physical  mesh  plane  distribution 
is  shown  in  Figure  3.  The  mesh  generation  has  been  applied  to  several  aft-swept 
configurations,  several  forward-swept  configurations  and  to  several  arbitrary 
wing-canard  parametric  variations.  Good  results  were  observed  in  all  cases.  For 
forward  swept  wings,  the  canard  tip  and  wing  tip  may  be  at  the  same  streamwise 
location.  In  this  case,  the  two  streamwise  points  used  for  £  =  ±  1  in  the 
transformation  are  the  wing/canard  tip  location  and  near  the  trailing  edge  of  the 
centerline  of  the  wing. 

The  spanwise  grid  transformation  is  much  simpler: 

Y  =  (l/A2)TANH"1(n)  +  CjH  >  C3n3  (8) 

This  transforms  the  domain  0<£<1  to  0<Y<°°.  The  Constant  A„  is  determined  to  place 
18  spanwise  mesh  planes  between  the  centerline  and  wing  tip,  with  the  wing  tip 
falling  midway  between  computational  planes.  The  constants  C1  and  C~  are  used  to 
perturb  the  grid  so  that  the  canard  tip  falls  midway  between  two  mesn  planes  and 
the  wing  tip  location  is  not  changed.  The  constants  and  are  set  to  zero  when 
a  canard  is  not  present. 

The  vertical  grid  transformation  is 

Z  =  A3TAN((n/2)t)  +  DjC  +  D3£3  (9) 

This  transforms  - 1 to  _®£Z<®.  The  constant  A,  is  chosen  to  provide  sufficient 
mesh  stretching  for  proper  supersonic  flow  development.  The  constants  and  D., 
are  set  to  make  the  canard  plane  and  a  mesh  plane  coincident,  while  maintaining  a 
monotonic,  positive  grid  stretching  (i.e.,  {h^n/2  +  Dj)  >  0).  When  a  canard  is  not 
present,  the  constants  and  D3  are  set  to  zero. 
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3.  MULTIPLE  WAKE  ANALYSIS 


With  two  lifting  surfaces  and  their  attendant  wakes,  the  basic  algorithm  for 
the  SLOR  solution  was  modified.  The  current  capability  of  the  code  is  an  above 
wing,  non-overlapping  canard.  As  sketched  in  Figure  2,  nine  distinct  flow-field 
regions  for  vertical  line  relaxation  occur: 

Region  1  -  Bounded  by  top  and  bottom  infinity 
Region  2  -  Bounded  by  the  canard  and  top  infinity 
Region  3  -  Bounded  by  bottom  infinity  and  the  canard 

Region  4  -  Bounded  by  the  canard  wake  and  top  infinity 

Region  5  -  Bounded  by  bottom  infinity  and  the  canard  wake 

Region  6  -  Bounded  by  the  wing  and  canard  wake 

Region  7  -  Bounded  by  bottom  infinity  and  the  wing 

Region  8  -  Bounded  by  the  wing  wake  and  the  canard  wake 
Region  9  -  Bounded  by  bottom  infinity  and  the  wing  wake 

A  wing  alone  solution  uses  only  vertical  line  regions  1  through  5.  The 
potential  jump  in  the  wake  is  constant,  equal  to  that  at  the  trailing  edge  of  the 
lifting  surface.  As  indicated  in  Figure  2,  the  wakes  are  constrained  to  remain  in 
the  plane  of  the  corresponding  lifting  surface.  This  is  considered  to  be 
consistent  with  the  other  small  disturbance  approximations  in  the  solution,  but 
actually  is  also  done  in  full  potential  codes.  When  the  wake  of  an  upstream 
lifting  surface  comes  arbitrarily  close  to  a  second  surface,  the  error  of  the 
undeflected  wake  model  may  be  significant.  Unfortunately,  detailed  pressure  data 
for  simple  wing-canard  geometries  does  not  exist  to  isolate  this  effect  and 
evaluate  the  numerical  model.  Good  force  and  wing  pressure  predictions  were 
obtained  for  the  limited  analyses  down  in  this  study. 

4.  EMBEDDED  FINE  GRID  SYSTEM 

Individual  fine  grid  arrays  are  constructed  for  the  wing  and  canard.  These 
secondary  mesh  systems  serve  two  purposes.  First,  detailed  computations  are 
performed  only  in  a  region  very  close  to  the  surface  where  gradients  are  large  and 
details  are  important.  The  resulting  numerical  efficiency  permits  a  very  dense 
computational  mesh,  a  benefit  in  both  the  resolution  of  shock  waves  and  the 
calculation  of  configuration  forces  and  moments.  Second,  the  embedded  mesh 
systems  are  independent  and  optimized  for  a  particular  geometric  component.  The 
system  is  not  constrained  by  a  single  geometry-fitting  transformation.  This  will 
facilitate  future  applications  to  configurations  with  multiple  components. 

Fine  grid  arrays  are  set  up  at  each  position  where  a  crude  spanwise  mesh 
plane  cuts  the  wing  and  canard  surface.  With  18  spanwise  crude  planes  between  the 
centerline  and  wing  tip,  an  isolated  wing  analysis  would  have  18  fine  wing  grids. 
(A  fuselage  would  reduce  the  number  of  fine  wing  grids  by  the  number  of  mesh  planes 
within  the  computational  wing  root  junction.)  The  number  of  canard  fine  grids  is 
proportional  to  the  extent  of  the  canard  semispan.  If  the  canard  semispan  were 
half  that  of  the  wing,  then  the  canard  would  have  approximately  half  the  number  of 
wing  fine  grids.  This  would  increase  the  computation  time  for  the  fine  grid 
solution  by  fifty  percent.  The  fine  grid  arrays  are  evenly  spaced  in  both  the 
streamwise  and  vertical  directions. 
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The  local  section  leading  edge  is  placed  midway  between  mesh  points  and  the 
trailing  edge  is  placed  at  a  mesh  point.  The  vertical  mesh  spacing  is  the  same  for 
all  wing  planes.  The  streamwise  mesh  spacing  is  scaled  by  the  local  chord.  Thus 
each  section  has  the  same  number  of  grid  points  along  the  chord.  Typically  80  mesh 
points  are  placed  on  the  chord,  but  code  dimensions  will  allow  up  to  100  points  on 
the  chord  if  greater  resolution  is  desired.  The  vertical  spacing  is  adjusted  to 
the  appropriate  value  for  the  canard  fine  grids.  As  for  the  wing,  the  vertical 
mesh  spacing  is  the  same  for  all  canard  planes,  and  the  streamwise  spacing  is 
scaled  to  the  local  chord.  The  boundaries  of  the  fine  grids  are  at  20%  of  the 
local  chord  in  front  of  each  leading  edge  and  10%  behind  each  trailing  edge.  The 
upper  boundary  is  at  30%  of  the  wing  average  chord  and  the  lower  boundary  at  10%. 
The  total  number  of  fine  grid  points  will  vary  for  different  canard/wing  semispan 
ratios.  The  dimensions  of  the  code  will  allow  108,000  fine  grid  points.  This 
corresponds  to  32  fine  grids  on  the  wing  and  canard,  with  each  having  135 
streamwise  points  and  25  vertical  points. 

The  fine  grid  system  is  illustrated  in  Figure  3.  As  described  above,  the 
fine  grids  are  sheared  and  tapered  to  conform  to  the  wing  and  canard  planforms. 
This  represents  a  transformation  function  between  the  fine  grid  computational 
domain  (4,n>0  an<^  the  physical  domain  (x,y,z): 

x  =  C(XGTE  -  XGLE)  +  XGLE 

tanh'^g) 

y  C 

L2 

z  =  £  (10c) 

The  spanwise  (y)  transformation  is  the  same  as  that  for  the  crude  grid.  The 
embedded  grid  approach  allows  independent  transformations  for  the  wing  and  canard 
in  the  streamwise  (x)  and  vertical  (z)  directions.  Thus  the  sweep  and  taper  of  the 
wing  does  not  effect  the  canard  grid  system  and  the  canard  sweep  and  taper  does  not 
effect  the  wing  grid  system. 

The  wing  fine  grid  solution  need  allow  for  the  presence  of  the  canard  wake. 
Figure  4  illustrates  the  vertical  line  relaxation  regions  for  the  wing  fine  grid 
that  includes  the  canard  wake.  The  different  regions  have  their  obvious  counter¬ 
part  in  the  crude  grid  solution  approach  (Figure  2).  Within  an  individual  fine 
grid,  the  vertical  line  relaxation  starts  at  the  upstream  boundary  and  moves  to 
the  downstream  boundary.  The  fine  grids  are  solved  sequentially,  beginning  with 
the  wing  root.  The  solution  proceeds  to  the  wing  tip  and  then  begins  the  canard 
solution  at  the  canard  root  juncture.  When  the  fine  grid  solution  at  the  canard 
tip  has  been  updated,  the  fine  grid  solution  iteration  is  completed. 

As  shown  in  Figure  4,  the  wing  fine  grid  solution  allows  only  the  canard 
wake  -  not  the  canard  surface  -  above  the  wing  surface.  The  analysis  is  restricted 
to  non-overlapping,  high  canard  planforms.  Strictly,  the  wing  (or  canard)  surface 
must  not  penetrate  the  canard  (or  wing)  fine  grid  boundary.  A  mild  violation  of 
this  restriction  does  not  disrupt  the  solution  development.  A  significant 
violation  will  lead  to  rapid  divergence  of  the  fine  grid  solution. 


(10a) 

(10b) 
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5.  FINE/CRUDE  GRID  INTERFACE 

The  embedded  fine  grid  interface  with  the  crude  grid  is  accomplished  by 
alternately  updating  the  crude  grid  and  fine  grid  solutions.  Potential  values  at 
interior  points  (i.e.,  on  a  lifting  surface)  of  the  crude  grid  are  fixed  by 
interpolating  the  most  recent  fine  grid  solution.  The  potential  values  at  the 
perimeter  of  the  fine  grids  are  fixed  by  interpolating  the  most  recent  crude  grid 
solution.  In  order  to  speed  up  the  overall  solution  convergence,  the  fine  grid 
solution  is  not  calculated  until  the  crude  grid  solution  has  established  the 
"coarse"  characteristics  of  the  flow  (typically  75  iterations).  Then  the 
crude/fine  interaction  is  begun  until  both  grids  are  satisfactorily  converged 
(typically  an  additional  95  cycles).  Linear  interpolation  of  the  crude  grid 
potential  values  is  used  to  initialize  the  fine  grid. 

One  cycle  of  the  crude/fine  grid  interaction  consists  of  two  steps  (see 
Figure  5) : 

1.  The  embedded  wing  grid  is  swept  holding  fine  grid  perimeter  points 
fixed  as  an  outer  boundary.  Conventional  Neumann  boundary 
conditions  (<|>  )  are  imposed  at  fine  grid  section  boundary  points 
forming  an  inner  boundary. 

2.  The  crude  grid  section  boundary  points  are  computed  using  the 
potentials  at  the  fine  grid  section  boundary  points  (linear 
interpolation).  These  crude  potentials  (<p)  are  held  fixed  for  the 
global  crude  grid  sweep  forming  an  array  of  Dirichlet  inner 
boundary  conditions.  Infinity  boundary  conditions  at  the  limits 
of  the  crude  computational  space  form  the  outer  boundary.  At  the 
end  of  the  crude  grid  sweep,  crude  grid  potentials  are  used  to 
update  the  fine  grid  perimeter  points. 

One  solution,. iteration  of  only  the  crude  grid  takes  1.4  seconds  of  CDC  7600 
CPU  time  (3.4x10  CPU  seconds  per  point).  One  crude/fine  cycle  takes  4.2  seconds 
CPU  time  for  23  wing  and  canard  fine  grids.  Each  fine  grid  has  120  streamwise 
points  and  25  vertical ^points .  Thus  one  crude/fine  solution  cycle  solves  110,000 
points  total  (3.8x10  CPU  seconds  per  point).  For  different  wing-canard 
combinations,  the  crude/fine  CPU  time  is  essentially  proportional  to  the  number  of 
wing  and  canard  fine  grid  arrays. 

The  development  of  the  solution  is  monitored  by  calculating  the  maximum 
update  to  the  flow  field  potential  (DPM  in  the  printout)  and  its  position  in  the 
three-dimensional  flow  field.  Typically,  for  most  configurations  at  flow  condi¬ 
tions  which  are  of  interest,  the  value  will  start  near  5  x  10  and  end  near  5  x 
10  for  the  final  solution.  The  solution  is  essentially  converged  when  the 
circulation  (CIR  in  the  printout)  shows  no  significant  change  in  about  twenty 
solution  cycles.  It  is  possible  that  a  non-decreasing  DPM  value  is  due  to  "noise" 
in  the  crude/fine  interface  that  does  not  decay.  A  two  or  three  order  of  magnitude 
DPM  reduction  should  be  expected.  The  position  of  DPM  is  useful  in  pinpointing 
the  problem  area  when  convergence  is  hindered.  Problem  areas  may  develop  when  the 
flow  conditions  are  extreme  or  when  the  geometric  representation  is  in  error. 
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Three  other  parameters  that  indicate  the  solution  convergence  appear  in  the 
code  output  as  DPMAV,  RSD  and  RSDAV.  The  DPMAV  parameter  is  the  average  of  the 
absolute  value  of  the  flow  potential  corrections  at  every  point.  It  is  a  more 
reliable  indication  of  the  solution  convergence  than  DPM,  and  is  used  to  satisfy 
the  convergence  test  input  option.  The  RSD  parameter  is  the  largest  "error"  in 
the  solution.  The  location  of  the  RSD  value  is  also  printed.  This  "error"  is  the 
right  hand  side  of  the  finite  difference  equation  being  solved.  It  is  calculated 
with  both  "old"  and  "new"  flow  potential  values.  The  RSDAV  parameter  is  the 
average  of  the  absolute  value  of  the  error  at  every  point. 

The  indices  for  the  location  of  DPM  and  RSD  are  shifted  when  the  point  is 
located  at  the  lower  side  of  the  wing  or  canard  plane.  The  jump  in  potential  at 
these  locations  is  handled  by  an  extra  array  of  potential  values  that  have  their 
own  index  system.  In  the  crude  grid,  500  is  added  to  the  streamwise  index  for  the 
canard  plane  lower  side  location  and  500  +  IMAX  is  added  to  the  streamwise  index 
for  the  wing  plane  lower  side  location  (the  default  for  IMAX  is  51).  In  the  fine 
grid  system,  500  is  added  to  the  streamwise  index  for  the  plane  of  the  lifting 
surface  (wing  in  a  wing  fine  grid,  canard  in  a  canard  fine  grid)  and  500  +  IMAXW  is 
added  to  the  streamwise  index  for  a  canard  wake  above  the  wing  surface  (the 
default  for  IMAXW  is  120). 

No  precise  non-dimensional  parameter  is  available  to  compare  the  solution 
convergence  for  different  grids.  The  DPM  and  DPMAV  values  are  scaled  by  the  wing 
average  cord.  This  results  in  the  same  values  for  DPM  and  DPMAV  if  the  scale  of 
the  input  geometry  is  changed  (e.g.,  full  scale  or  model  scale).  Changes  to  the 
planform  (e.g.,  aspect  ratio  or  taper)  would  result  in  DPM  and  DPMAV  values  that 
are  not  directly  comparable.  As  mentioned  above,  the  reduction  in  DPM  and  DPMAV 
relative  to  the  starting  value  should  be  used  to  evaluate  the  solution 
convergence. 

6.  BOUNDARY  LAYER  CALCULATION 

Viscous  effects  are  computed  in  the  analysis  code  by  coupling  a  modified 
Bradshaw  boundary  layer  computation  with  the  inviscid  potential  flow  solution. 
The  boundary  layer  calculation  is  virtually  identical  to  the  method  developed  by 
Mason  [2],  The  method  employs  the  modified  chord  technique  of  Nash  [3],  which 
represents  an  infinite  sheared  wing  boundary  layer  calculation.  The  wing  sweep 
angle  is  that  of  the  local  mid-chord  span  line,  such  that  it  may  vary  across  the 
span.  The  2-D  Bradshaw  turbulent  boundary  layer  analysis  [4]  provides  the 
foundation  for  the  method.  The  use  of  a  modified  2-D  boundary  layer  analysis 
greatly  reduces  the  necessary  computer  time  and  has  demonstrated  good  results  for 
several  different  codes  [1,2,5]. 

The  boundary  layer  calculation  provides  a  displacement  thickness  and  skin 
friction  calculation  at  each  analysis  station  of  the  wing  and  canard.  The  slope 
of  the  displacement  thickness  is  used  to  modify  the  surface  boundary  conditions  in 
the  inviscid  solution.  The  local  skin  friction  calculation  is  used  to  provide  a 
viscous  drag  estimate  for  the  configuration  at  the  end  of  the  analysis  run. 
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Surface  pressures  are  calculated  from  the  fine  grid  potential  solution  every 
20  cycles.  The  boundary  layer  calculation  for  twenty  three  wing-canard  stations 
(an  upper  and  lower  surface  calculation  at  each  station)  takes  approximately  22 
seconds  of  CDC  7600  CPU  time.  Thus  four  boundary  layer  calculations  during  the 
inviscid  solution  with  a  boundary  layer  calculation  at  the  end  will  add  110 
seconds  of  CPU  time  to  that  for  the  inviscid  solution  alone.  As  mentioned  for  the 
inviscid  solution,  the  viscous/inviscid  interaction  is  essentially  converged  when 
the  change  in  lift  (circulation)  is  negligible. 

7.  FUSELAGE  MODELING 

The  body  is  modeled  in  the  solution  by  a  constant  cross-section  computational 
surface  in  the  Cartesian  crude  grid.  The  input  data  allow  for  simple  axisymmetric 
body  definition  or  detailed  "Quick  Geometry"  [6,7]  body  definition.  Body  boundary 
conditions  are  imposed  by  fixing  the  velocity  potential  values  on  the  body 
computational  surface.  The  procedure  follows  that  of  Reference  [2]  and  is 
described  in  detail  in  Reference  [1].  Body  pressures  from  equation  (3)  are  used 
to  produce  a  calculation  of  the  body  force  and  moment  contribution. 

SECTION  III 

COPES/CONMIN  OPTIMIZATION 

The  wing-body-canard  analysis  code  was  coupled  with  the  COPES  and  CONMIN 
routines  of  Vanderplaats  [8,9].  The  COPES  code  is  a  control  program  that  connects 
the  numerical  optimization  code  CONMIN  with  the  aerodynamic  analysis  code.  The 
structure  of  the  PANDORA  code  allows  the  numerical  optimization  to  be  coupled  to 
any  analysis  code.  Several  "analysis"  codes  could  be  coupled  together  to  provide 
information  for  the  optimization.  This  information  might  be  supersonic 
performance,  TOGW  changes  or  mechanical  system  sizing  estimates.  The  available 
computer  resources  represent  the  only  limit  of  complexity. 

The  COPES  and  CONMIN  routines  were  slightly  modified  for  inclusion  in  the 
PANDORA  code.  Changes  were  made  in  the  main  COPES  routine  (Program  COPES  in 
0verlay(0 ,0))  and  subroutine  CNMN06  (in  0verlay(5 ,0)) .  The  change  to  CNMN06  is 
important.  It  provides  subroutine  ANALIZ  with  the  index  of  the  best  result  during 
an  optimization  search.  In  this  way,  the  next  search  solution  can  be  restarted 
from  the  best  previous  result.  This  more  closely  models  the  search  strategy  in 
CONMIN,  resulting  in  more  accurate  information  from  the  aerodynamic  analysis.  The 
changes  in  the  COPES  and  CONMIN  routines  are  identified  in  the  code  with  comment 
cards . 

The  computer  code  is  written  in  FORTRAN,  employing  the  CDC  overlay  structure. 
Storage  requirements  on  a  CDC  7600  are  140K„  small  core  memory  and  770Kg  large 
core  memory.  Typical  CPU  times  on  the  7600  for  a  wing-body-canard  analysis  with 
viscous  effects  is  eleven  minutes.  The  CPU  time  for  an  optimization  run  will  vary 
according  to  the  complexity  of  the  optimization  problem,  as  discussed  below. 

The  optimization  algorithm  within  CONMIN  is  a  modified  Method  of  Feasible 
Directions  [9].  The  gradient  information  for  the  algorithm  is  calculated  by 
sequentially  perturbing  each  design  variable.  Each  design  variable  perturbation 
requires  analysis  by  the  flow  solution  routines.  Thus,  the  computer  time  is 


8 


proportional  to  the  number  of  design  variables.  The  gradient  information  estab¬ 
lishes  a  search  direction  that  should  improve  the  design  (decrease  the  objective 
function)  while  satisfying  any  constraints.  The  search  direction  is  "explored" 
until  a  relative  optimum  is  found  or  any  constraints  become  violated.  If  the 
starting  conditions  violated  any  constraints,  the  search  direction  will  be  that 
which  satisfies  the  constraints  with  the  least  objective  function  increase.  One 
to  four  flow  solutions  are  required  during  the  search.  Completion  of  the  search 
constitutes  one  optimization  iteration.  The  following  discussion  of  COPES  usage 
is  excerpted  from  Reference  [8].  A  more  detailed  description  will  be  found  in  the 
reference.  Reference  [9]  should  be  consulted  if  more  details  of  the  numerical 
optimization  algorithm  is  desired. 

If  it  is  desired  to  run  only  a  simple  analysis  using  COPES,  only  three  data 
cards  are  required  for  the  COPES  program:  a  TITLE  card,  a  control  parameter, 
(NCALC  =  1),  and  an  END  card.  If  the  optimization  or  parametric  analysis  (sensi¬ 
tivity)  capabilities  of  COPES  are  to  be  used,  additional  data  must  be  read.  This 
data  will  identify  which  parameters  in  the  global  common  block,  GLOBCM,  are  used. 
To  set  up  the  COPES  data,  the  user  must  have  a  basic  understanding  of  how  the  data 
in  the  global  common  block  is  accessed  by  COPES.  This  is  outlined  in  the  following 
section. 

1.  OPTIMIZATION  DATA  MANAGEMENT 

In  order  to  perform  design  operations,  the  COPES  program  must  access  the  data 
in  common  block  GLOBCM.  This  done  by  defining  the  location  in  GLOBCM  where  a 
specified  parameter  resides.  For  example,  consider  the  common  block  for  a 
cantilevered  beam  design  problem: 

COMMON/ GLOBCM/ B , H , VOL , BSTRES , SHRSTR , DELTA , HB , E , AL 

The  volume  of  material,  VOL,  is  the  third  parameter  in  the  common  block;  that 
is,  it  resides  in  location  3,  referred  to  as  the  global  location  number. 
Similarly  the  bending  stress,  BSTRES,  is  in  global  location  4  and  the  beam  width 
is  in  global  location  1.  Thus,  the  parameters  are  referred  to  by  their  respective 
location  numbers  in  global  common. 

For  convenience  in  preparing  data  for  the  COPES  program,  a  simple  "CATALOG" 
of  parameters  may  be  defined.  For  the  cantilevered  beam,  this  catalog  would  be: 

GLOBAL  FORTRAN 

LOCATION  NAME  DEFINITION _ 


1 

B 

Beam  width 

2 

H 

Beam  height 

3 

VOL 

Volume  of  material 

4 

BSTRES 

Maximum  bending  stress 

5 

SHRSTR 

Maximum  shear  stress 

6 

DELTA 

Deflection  under  the  load 

7 

HB 

Ratio,  H/B 

8 

E 

Young's  modulus 

9 

AL 

Length  of  beam 

9 


As  another  example,  consider  a  global  common  block  containing  arrays: 


GLOBAL 

LOCATION 

FORTRAN 

NAME 

DEFINITION 

1 

A 

Area 

2 

Y(10) 

Vector 

12 

Q 

. 

13 

C(2,2) 

. 

17 

H 

etc. 

The  dimensions  are  given  with  the  FORTRAN  name  as  a  reminder  that  the 
parameter  is  an  array.  In  this  case,  the  third  parameter  in  the  Y  array  is  in 
global  location  4.  Remembering  that  arrays  are  stored  column  by  column  the  C ( 1 , 2 ) 
array  location  is  in  global  location  15. 

It  will  be  seen  that  identifying  parameters  according  to  their  location  in 
GLOBCM  provides  a  great  deal  of  flexibility  in  using  the  COPES  program  for  design. 

In  the  following  section,  definitions  of  terms  commonly  used  in  automated 
design  are  given  for  easy  reference. 

2.  OPTIMIZATION  TERMINOLOGY 

The  COPES  program  currently  provides  six  specific  capabilities: 

1.  Simple  analysis,  just  as  if  COPES  was  not  used. 

2.  Optimization  -  Minimization  or  maximization  of  one  calculated 
function  with  limits  imposed  on  other  functions. 

3.  Sensitivity  analysis  -  the  effect  of  changing  one  or  more  design 
variables  on  one  or  more  calculated  functions. 

4.  Two-variable  functions  space  -  analysis  for  all  specified  combi¬ 
nations  of  two  design  variables. 

5.  Optimum  sensitivity  -  same  as  sensitivity  analysis  except,  at  each 
step,  the  design  is  optimized  with  respect  to  the  remaining 
independent  design  variables. 

6.  Approximate  optimization  -  optimization  using  approximation 
techniques.  Usually  more  efficient  than  standard  optimization  for 
up  to  10  design  variables  of  if  multiple  optimizations  are  to  be 
performed. 

In  defining  the  data  required  to  execute  the  COPES  program,  the  following 
definitions  are  useful. 


Design  Variables  -  Those  parameters  which  the  optimization  program  is 
allowed  to  change  in  order  to  improve  the  design.  Design  variables  appear  only  on 
the  right  hand  side  of  equations  in  the  analysis  program.  COPES  considers  two 
types  of  design  variables,  independent  and  dependent.  If  two  or  more  variables 
are  always  required  to  have  the  same  value  or  be  in  a  constant  ratio,  one  is  the 
independent  variable  while  the  remaining  are  dependent  variables.  For  example,  if 
the  height  is  required  to  be  10  times  the  width  of  the  cantilevered  beam,  B  would 
be  the  independent  variable  while  H  would  be  the  dependent  variable. 

Objective  Function  -  The  parameter  which  is  to  be  minimized  or  maximized 
during  optimization.  Also  the  parameters  calculated  as  functions  of  specified 
design  variables  during  a  sensitivity  or  two-variable  function  space  study. 
Objective  functions  always  occur  on  the  left  side  of  equations,  unless  the 
objective  function  is  also  a  design  variable  (the  beam  height  may  be  minimized  as 
an  objective  function  if  it  is  also  a  design  variable.  In  this  way,  the  minimum 
height  is  found  for  which  no  constraints  are  violated).  An  objective  function  may 
be  linear  or  non-linear,  implicit  or  explicit,  but  must  be  a  function  of  the 
design  variables  to  be  meaningful. 

Constraint  -  Any  parameter  which  must  not  exceed  specified  bounds  for  the 
design  to  be  acceptable.  Constraint  functions  always  appear  on  the  left  side  of 
equations.  Just  as  for  objective  functions,  constraints  may  be  linear  or  non¬ 
linear,  implicit  or  explicit,  but  must  be  functions  of  the  design  variables. 

Constraint  Set  -  A  group  of  constraints  which  appear  consecutively  in  the 
global  common  block  and  which  all  have  the  same  limits  imposed.  This  is  a 
convenience  which  allows  several  constraints  to  be  identified  with  a  minimum  of 
data. 


Global  Common  -  Common  block  GLOBCM  containing  design  information. 

Global  Location  -  Location  of  a  particular  parameter  in  GLOBCM. 

3.  OPTIMIZATION  INPUT  DATA  FORMAT 

In  order  to  execute  the  COPES  program  it  is  necessary  to  provide  formatted 
data  for  COPES,  followed  by  data  for  the  ANALIZ  program  which  is  coupled  to  COPES. 
Section  4.1  defines  the  data  which  is  required  by  COPES.  The  data  is  segmented 
into  "BLOCKS"  for  convenience.  All  formats  are  alphanumeric  for  TITLE,  and  END 
cards,  F10  for  real  data  and  110  for  integer  data. 

The  COPES  data  begins  with  a  TITLE  card  and  ends  with  an  END  card.  This  is 
followed  by  data  to  be  read  by  the  user  supplied  subroutine  ANALIZ. 

Comment  cards  may  be  inserted  anywhere  in  the  COPES  data  stack  prior  to  the 
END  card,  and  are  identified  by  a  dollar  sign  ($)  in  column  1. 
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While  the  input  description  defines  COPES  data  in  formatted  fields  of  ten, 
the  data  may  actually  be  read  in  more  conveniently  by  separating  data  by  commas  or 
one  or  more  blanks.  If  more  than  one  number  is  contained  on  an  unformatted  data 
card,  a  comma  must  appear  somewhere  on  the  card.  If  exponential  numbers  such  as 
2.5+10  are  read  on  an  unformatted  card,  there  must  be  no  embedded  blanks . 
Unformatted  cards  may  be  intermingled  with  formatted  cards.  Real  numbers  on  an 
unformatted  card  must  have  a  decimal  point. 

Examples 

Unformatted  data: 

5, 7, 1.3, 1.0+20,0, -5.1 
5,7,1.3,1.0+20, ,-5.1 
5  7  1.3  1.0+20,,  -  5.1 
5  7  1.3,  1.0+20  0  -5.1 

Equivalent  formatted  data: 

col  10  20  30  40  50  60  70  80 

5  7  1.3  1.0+20  0  -5.1 

Unformatted  data 

2 

2,3 
2  3 

Equivalent  formatted  data: 

col  10  20  30 

2 

2  3 

2  3  Note:  These  data  contain  no  commas,  so  it 

is  assumed  to  be  formatted  already. 


Unformatted  data: 

1.2.3.4.5.6.7.8.9.10.11 

col  10  20  30  40  50  60  70  80 

1  2  3  4  5  6  7  8 

9  10  11 

Note  that  two  formatted  data  cards  are  created  here.  Unformatted  data: 

1,2, 3, 4, 5, 6 

7.8.9.10.11 


12 


Equivalent  formatted  data: 

col _ 10 _ 20 _ 30 _ 40 _ 50 _ 60 _ 70 _ 80 

1  2  3  4  5  6 

7  8  9  10  11 

Note  that  the  above  two  examples  do  not  produce  the  same  formatted  data  cards. 

SECTION  IV 

INPUT  DATA  DESCRIPTION 

This  section  describes  the  necessary  solution/optimization  parameters  and 
geometry  input  needed  to  operate  the  PANDORA  code.  The  information  of  Section  III 
is  important  for  using  the  optimization.  Only  the  one-cycle  analysis  (NCALC  =  1 
in  DATA  BLOCK  B)  and  optimization  (NCALC  =  2)  COPES  options  have  been  used  in  this 
effort.  The  transonic  analysis  input  is  an  extension  of  that  from  the  base  code 
[l].  A  set  of  Usage  Notes  appears  at  the  end  of  the  analysis  input  description  to 
elaborate  on  the  less  familiar  input  parameters.  The  reader  is  reminded  that 
Volume  2,  Part  2  of  this  final  report  contains  important  discussions  of  the  code 
application  and  accuracy. 

1.  OPTIMIZING  INPUT  DATA 

DATA  BLOCP  A 

DESCRIPTION  Title  card. 

FORMAT  AND  EXAMPLE 

1  2  3  4  5  6  7  8  FORMAT 

TITLE  20A4 

CANTILEVERED  BEAM  DESIGN 

FIELD  CONTENTS 

1-8  Any  80  character  title  may  be  given  on  this  card. 
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DATA  BLOCK 


B 


DESCRIPTION: 

FORMAT  AND  EXAMPLE 

1  2 
NCALC  NDV 

2 

FIELD 

1  NCALC: 

0  - 

1  - 

2  - 

3  - 

4  - 

5  - 

6  - 

2  NDV: 

4  N2VAR: 

5  NXAPRX: 

6  IPNPUT: 

0  - 
1  - 
2  - 

7  IPDBG: 

0  - 
1  - 


Program  Control  Parameters. 


3 

4 

5 

6 

7  FORMAT 

NSV 

N2VAR 

NXAPRX  IPNPUT 

IPDBG 

7110 

2 

3 

5  2 

CONTENTS 

0 

0 

Calculation  Control 

Read  input  and  stop.  Data  of  blocks  A,  B  and  V  is  required. 
Remaining  data  is  optional. 

One  cycle  through  program.  The  same  as  executing  ANALIZ 
stand-along  (i.e.,  no  optimization).  Data  of  blocks  A,  B  and 
V  is  required.  Remaining  data  is  optional. 

Optimization.  Data  of  blocks  A-I  and  V  is  required.  Remain¬ 
ing  data  is  optional. 

Sensitivity  analysis.  Data  of  blocks  A,  B,  P,  Q  and  V  is 
required.  Remaining  data  is  optional. 

Two  variable  function  space.  Data  of  blocks  A,  B,  and  R-V  is 
required.  Remaining  data  is  optional. 

Optimum  Sensitivity.  Data  of  blocks  A-I,  P,  Q,  and  V  is 
required.  Remaining  data  is  optional. 

Optimization  using  approximation  techniques.  Data  of  blocks 
A-0  and  V  is  required.  Remaining  data  is  optional 

Number  of  design  variables  on  which  sensitivity  analysis  or 
optimization  will  be  performed. 

Number  of  objective  functions  in  a  two  variable  function 
space  study. 

Number  of  X-variables  for  approximate  analysis/optimization. 
>  NDV 

Input  print  control 

Print  card  images  of  data  plus  formatted  print  of  input  data 
Formatted  print  only  of  input  data 
No  print  of  input  data. 

Debug  print  control 
Off 

On,  ANALIZ  called  for  output  after  each  analysis. 
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DATA  BLOCK 


C  OMIT  IF  NDV  =  0  IN  BLOCK  B 


DESCRIPTION:  Integer  optimization  control  parameter. 

FORMAT  AND  EXAMPLE 


1 

2 

3 

4 

5  6 

7 

8 

FORMAT 

IPRINT 

ITMAX 

ICNDIR 

NSCAL 

ITRM 

LINOBJ 

NACMX1 

NFDG 

5  0 

0 

0 

0 

0 

0 

0 

FIELD 


CONTENTS 


1 


IPRINT:  Print  control  used  in  the  optimization  program  CONMIN. 

0  -  No  print  during  optimization. 

1  -  Print  initial  and  final  optimization  information. 

2  -  Print  above  plus  objective  function  value  and  design  variable 

values  at  each  iteration. 

3  -  Print  above  plus  constraint  values,  direction  vector  and  move 

parameter  at  each  iteration. 

4  ~  Print  above  plus  gradient  information. 

5  -  Print  above  plus  each  proposed  design  vector,  objective 

function  and  constraint  values  during  the  one-dimensional 
search. 


2 

3 

4 


5 


6 


7 


ITMAX:  Maximum  number  of  optimization  iterations  allowed.  DEFAULT  = 

20. 

ICNDIR:  Conjugate  direction  restart  parameter.  DEFAULT  =  NDV  +  1. 

NSCAL :  Scaling  parameter.  .GT.O  -  Scale  design  variables  to  order  of 

magnitudes  one  every  NSCAL  iterations.  .LT.O  -  Scale  design 
variables  according  to  user-input  scaling  values.  Good  value 
is  ICNDIR,  1  is  not. 

ITRM:  Number  of  consecutive  iterations  which  must  satisfy  relative 

or  absolute  convergence  criterion  before  optimization  process 
is  terminated.  DEFAULT  =  3. 

LINOBJ:  Linear  objective  function  identifier.  If  the  optimization 

objective  is  known  to  be  a  linear  function  of  the  design 
variables,  set  LINOBJ  =  1.  DEFAULT  =  Non-linear. 

NACMX1:  One  plus  the  maximum  number  of  active  constraints  antici¬ 

pated.  DEFAULT  =  NDV  +  2. 
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FIELD  (coat'd) 


CONTENTS  (coat'd) 


8 


NFDG:  Fiaite  difference  gradient  identifier. 

0  -  All  gradient  information  is  computer  by  finite  difference 

with  CONMIN. 

1  -  All  gradient  information  is  computed  analytically  by  the 

user-supplier  code. 

2  -  Gradient  of  objective  is  computed  analytically.  Gradients  of 

constraints  are  computed  by  finite  difference  within  CONMIN. 


REMARKS 

1)  Currently  NFDG  must  be  zero  in  COPES. 

2)  IPRINT  =  5  is  recommended. 
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DATA  BLOCK 


D  OMIT  IF  NDV  =  0  IN  BLOCK  B 


DESCRIPTION :  Floating  point  optimization  program  parameters. 

FORMAT  AND  EXAMPLE 

1  2  3  4  5  6  7  FORMAT 

FDCH  FDCHM  CT  CTMIN  CTL  CTLMIN  THETA  7F10 

0.0  0.0  0.0  0.0  0.0  0.0  0.0 

DELFUN  DABFUN  ALPHAX  AB0BJ1 

0.0  0.0  0.0  0.0 

NOTE:  TWO  CARDS  ARE  READ  HERE . 

FIELD  CONTENTS 

1  FDCH:  Relative  change  in  design  variables  in  calculating  finite 

difference  gradients.  DEFAULT  =0.01 

2  FDCHM:  Minimum  absolute  step  in  finite  difference  gradient  calcula¬ 

tions.  DEFAULT  =  0.001. 

3  CT:  Constraint  thick  parameter.  DEFAULT  =0.05.  s 

4  CTMIN:  Minimum  absolute  value  of  CT  considered  in  the  optimization 

process.  DEFAULT  =  0.004. 

5  CTL:  Constraint  thickness  parameter  for  linear  constraints. 

DEFAULT  =  -0.01. 

6  CTLMIN:  Minimum  absolute  value  of  CTL  considered  in  the  optimization 

process.  DEFAULT  =  0.001. 

7  THETA:  Mean  value  of  push-off  factor  in  the  method  of  feasible 

directions.  DEFAULT  =1.0. 

1  DELFUN:  Minimum  relative  change  in  objective  function  to  indicate 

convergence  of  the  optimization  process.  DEFAULT  =  0.001. 

2  DABFUN:  Minimum  absolute  change  in  objective  function  to  indicate 

convergence  of  the  optimization  process.  DEFAULT  =  0.001 
times  the  initial  objective  value. 

3  ALPHAX:  Maximum  fractional  change  in  any  design  variable  for  first 

estimate  of  the  step  in  the  one-dimensonal  search.  DEFAULT  = 

0.1. 
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FIELD  (cont'd) 

4  AB0BJ1 

REMARKS 

T5  The  DEFAULT 


CONTENTS  (cont'd) 

:  Expected  fractional  change  in  the  objective  function  for 

first  estimate  of  the  step  in  the  one-dimensional  search. 
DEFAULT  =0.1 

values  for  these  parameters  usually  work  well. 


18 


DATA  BLOCK 

E 

OMIT  IF 

NDV 

=  0  IN  BLOCK  B 

DESCRIPTION 

:  Total  number 

of  design  variables, 

design  objective 

identification 

and 

sign. 

FORMAT  AND 

EXAMPLE 

1 

2 

3 

FORMAT 

NDVTOT 

IOBJ 

SGNOPT 

2110, F10 

0 

3 

-1.0 

FIELD  CONTENTS 


1  NDVTOT:  Total  number  of  variables  linked  to  the  design  variables. 

This  option  allows  two  or  more  parameters  to  be  assigned  to  a 
single  design  variable.  The  value  of  each  parameter  is  the 
value  of  the  design  variable  times  a  multiplier,  which  may  be 
different  for  each  parameter.  DEFAULT  =  NDV. 

2  IOBJ:  Global  variable  location  associated  with  the  objective 

function  in  optimization. 

3  SGNOPT:  Sign  used  to  identify  whether  function  is  to  be  maximized  or 

minimized.  +1.0  indicates  maximization.  -1.0  indicates 
minimization.  If  SGNOPT  is  not  unity  in  magnitude  it  scales 
the  magnitude  of  the  objective. 
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DATA  BLOCK 


OMIT  IF  NDV  =  0  IN  BLOCK  B 


DESCRIPTION: 


Design  variable  bounds,  initial  values  and  scaling  factors. 


FORMAT  AND  EXAMPLE 


FORMAT 

4F10 


NOTE:  READ  ONE  CARD  FOR  EACH  OF  THE  NDV  INDEPENDENT  DESIGN  VARIABLES.  VALUES 

ARE  IN  ORDER  FOR  DESIGN  VARIABLES  NDSGN  SEQUENCE  (I.E. ,  INPUT  SEQUENCE) 

FIELD  CONTENTS 

1  VLB:  Lower  bound  on  the  design  variable.  If  VLB.LT. -1 .OE+15 ,  no 

lower  bound. 

2  VUB:  Upper  bound  on  the  design  variable.  If  VUB.GT. 10.E+15 ,  no 

upper  bound. 

3  X:  Initial  value  of  the  design  variable.  If  X  is  non-aero,  this 

will  supercede  the  value  initialized  by  the  user-supplied 
subroutine  ANALIZ. 


SCAL: 


Design  variable  scale  factor.  Not  used  if  NSCAL.GE.O  in  BLOCK 


DATA  BLOCK 


G  OMIT  IF  NDV  =  0  IN  BLOCK  B 


DESCRIPTION:  Design  variable  identification. 

FORMAT  AND  EXAMPLE 


1 

2 

3 

NDSGN 

IDSGN 

AMULT 

1 

1 

1.0 

FORMAT 

2I10.F10 


NOTE:  READ  ONE  CARD  FOR  EACH  OF  THE  NDVTOT  DESIGN  VARIABLES. 


FIELD 


CONTENTS 


1  NDSGN:  Design  variable  number  associated  with  this  variable. 

2  IDSGN:  Global  variable  number  associated  with  this  variable. 


3 


AMULT 


Constant  multiplier  on  this  variable.  The  value  of  the 
variable  will  be  the  value  of  the  design  variable,  NDSGN, 
times  AMULT.  DEFAULT  =  1.0 
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DATA  BLOCK 


H  OMIT  IF  NDV  =  0  IN  BLOCK  B 


DESCRIPTION:  Number  of  constrained  parameters. 

FORMAT  AND  EXAMPLE 


1  FORMAT 

NCONS  110 

4 

FIELD 

CONTENTS 

1 

NCONS: 

Number  of  constraint  sets  in  the  optimization  problem 

REMARKS 


1)  If  two  or  more  adjacent  parameters  in  the  global  common  block  have  the  same 
limits  imposed,  these  are  part  of  the  same  constraint  set. 
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DATA  BLOCK 


I  OMIT  IF  NDV  =  0  IN  BLOCK  B.  OR  NCONS  =  0  IN  BLOCK  H 


DESCRIPTION:  Constraint  identification  and  constraint  bounds. 

FORMAT  AND  EXAMPLE 


1 

2 

3 

4 

FORMAT 

ICON 

JCON 

LCON 

3110 

4 

0 

0 

BL 

SCAL1 

BU 

SCAL2 

-1.0+20 

0.0 

20000. 

0.0 

NOTE:  READ  TWO  CARDS  FOR  EACH  OF  THE  NCONS  CONSTRAINT  SETS. 


FIELD  CONTENTS 


1 

ICON: 

First  global  number  corresponding  to  the  constraint  set. 

2 

JCON: 

Last  global  number  corresponding  to  the  constraint  set. 
DEFAULT  =  ICON. 

3 

LCON: 

Linear  constraint  identifier  for  this  constraint  set.  LCON 

1  indicates  linear  constraints.  If  in  doubt,  use  non-linear. 

1 

BL: 

Lower  bound  on  the  constrained  variables.  If  BL.LT. 
-1.0E+15,  no  lower  bound. 

2 

SCAL1: 

Normalization  factor  on  lower  bound.  DEFAULT  =  MAX  of 
(ABS(BL) ,  0.1). 

3 

BU: 

Upper  bound  on  the  constrained  variables.  If  BU.GT. 1 .0E+15 
no  upper  bound. 

4 

SCAL2: 

Normalization  factor  on  upper  bound.  DEFAULT  =  MAX  of 
(ABS(BU),  0.1). 

REMARKS 


1)  The  normalization  factor  should  usually  be  defaulted. 

2)  The  constraint  functions  sent  to  CONMIN  are  of  the  form:  (BL  -  VALUE)/  SCAL1 
.LE.  0.0  and  (VALUE  -  BU)/SCAL2  .LE.  0.0. 

3)  Each  constrained  parameter  is  converted  to  two  constraints  in  CONMIN  unless 
ABS(BL)  or  ABS(BU)  exceeds  1.0E+15,  in  which  case  no  constraint  is  created 
for  that  bound. 
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DATA  BLOCK  J  OMIT  IF  NXAPRX  =  0  IN  BLOCK  B 

DESCRIPTION:  Approximate  analysis/optimization  control  parameters. 

FORMAT  AND  EXAMPLE 

1  2  3  4  5  6  7  8  FORMAT 

NF  NXS  NXFS  NXA  INOM  ISCRX  IXCRXF  IPAPRX  8110 
5  1  1  1  0  0  0  1 

KMIN  KMAX  NPMAX  JNOM  INXLOC  INFLOC  6110 

0  0  0  0  0  0 


FIELD  CONTENTS 


1 

NF: 

Number  of  functions  to  be  approximated.  DEFAULT  =  number  of 
optimization  objective  and  constraint  functions. 

2 

NXS: 

Number  of  X-vectors  read  as  data. 

3 

NXFS: 

Number  of  X-F  pairs  read  as  data. 

4 

NXA: 

If  non-zero,  the  design  variables  read  by  SUBROUTINE  ANALIZ 
form  an  X-vector. 

5 

INOM: 

Nominal  X-vector  about  which  to  do  Taylor  expansion.  DEFAULT 
=  best  available. 

6 

ISCRX: 

File  from  which  NXS  X-vectors  are  read.  DEFAULT  =  5. 

7 

IXCRXF: 

File  from  which  NXFS  X-F  pairs  of  data  are  read.  DEFAULT  =  5. 

8 

IPAPRX: 

Print  control,  1  to  4.  4  is  most. 

1 

KMIN: 

Minimum  number  of  approximation  iterations. 

2 

KMAX: 

Maximum  number  of  approximation  iterations. 

3 

NPMAX: 

Maximum  number  of  designs  retained  for  Taylor  series 
expansion. 

4 


JNOM: 


Number  of  iterations  after  which  the  best  design  is  picked  as 
nominal . 


FIELD  (cont'd) 

5  INXLOC : 

6  INFLOC: 


CONTENTS  (cont'd) 

X-variable  global  location  identifier.  If  INXLOC  =  0,  the 
Taylor  series  expansion  is  on  the  design  variables  listed  in 
BLOCK  G.  , 

Function  global  location  identifier.  If  INFLOC  =  0,  the 
Objective  and  constraint  functions  identified  in  BLOCKS  E  and 
I  are  the  functions  on  which  the  Taylor  series  expansion  is 
performed. 


REMARKS 

1)  If  ISCRX  and/or  ISCRXF  file  number  is  other  than  5,  the  data  read  from  that 
file  is  assumed  to  be  binary  data. 

2)  If  NXS  =  NSFS  =  0,  NXA  is  defaulted  to  NXA  =  1,  even  if  it  is  read  as  zero. 
Also,  a  second  vector  of  design  variables  is  automatically  defined  by  COPES 
to  yield  two  independent  designs  to  start  the  optimization. 
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DATA  BLOCK 


K  OMIT  IF  NDV  =  0  IN  BLOCK  B,  OR  NXAPRX  =  0  IN  BLOCK  B 


DESCRIPTION :  Bounds  and  multipliers  for  approximate  optimization. 

FORMAT  AND  EXAMPLE 

1  2  3  4  5  6  7  8  FORMAT 

DX1  DX2  DX3  DX4  DX5  .  8F10 

.5  2. 

XFACT1  XFACT2  2F10 

0.  0. 


NOTE:  TWO  OR  MORE  CARDS  ARE  READ  HERE 


FIELD  CONTENTS 

1-8  DXI:  Allowable  change  (in  magnitude)  of  the  Ith  design  variable 

during  each  approximate  optimization. 

1  XFACT1:  Multiplier  on  DXI  when  the  diagonal  elements  of  the  H  matrix 

are  available.  DEFAULT  =1.5. 

2  XFACT2:  Multiplier  on  DFXI  when  all  elements  of  the  H  matrix  are 

available.  DEFAULT  =  2.0. 
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DATA  BLOCK 


OMIT  IF  NXAPRX  =  0  IN  BLOCK  B  OR  INXLOC  =  0  IN  BLOCK  J 


1-8  LOCI:  Global  location  of  Ith  ap.'-  mating  variable. 


REMARKS 

1)  If  INXLOC  =  0  in  BLOCK  J,  this  data  is  not  read.  In  this  case,  the  data  is 
defaulted  to  be  the  global  locations  of  the  design  variables  (IDSGN  values  in 
BLOCK  G) . 
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DATA  BLOCK  M  OMIT  IF  NXAPRX  =  0  IN  BLOCK  B  OR  INFLOC  =  0  IN  BLOCK  J 
DESCRIPTION :  Global  locations  of  functions  to  be  approximated. 

FORMAT  AND  EXAMPLE 

1  2  3  4  5  6  7  8  FORMAT 

LOCF1  LOCF2  LOCF3  LOCF4  .  8110 

3  5  6  4 


NOTE:  MORE  THAN  ONE  CARD  MAY  BE  READ  HERE. 


FIELD  CONTENTS 

1-8  LOCI:  Global  location  of  Ith  function  to  be  approximated. 

REMARKS 

1)  If  INFLOX  =  0  in  BLOCK  J,  this  data  is  not  read.  In  this  case,  the  data  is 
defaulted  to  be  the  global  locations  of  the  objective  function  (IOBJ  in  BLOCK 
E)  followed  by  the  global  locations  of  the  constrained  parameters  (ICON,  JCON 
in  BLOCK  I). 
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DATA  BLOCK 


N  OMIT  IF  NXS  =  0  IN  BLOCK  J 


DESCRIPTION:  X-vectors  for  approximate  optimization. 

FORMAT  AND  EXAMPLE 

12345678  FORMAT 

XII  XI 2  XI3  XI4  .  8F10 

4.  15. 

NOTE:  NXS  SETS  OF  DATA  ARE  READ  HERE. 

NOTE:  MORE  THAN  ONE  CARD  MAY  BE  READ  FOR  EACH  SET  OF  DATA. 

FIELD  CONTENTS 

1-8  XIJ:  Jth  value  of  Ith  X-vector,  J  =  1,NXAPRX. 
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DESCRIPTION:  X-F  pairs  of  information  for  approximate  optimization. 

FORMAT  AND  EXAMPLE 


2 

7200 


1 

2 

3 

4  5  6  7 

8 

FORMAT 

XI 

X2 

18. 

X3 

X4  . 

•  • 

8F10 

Yl 

i. 

Y2 

416.667 

Y3 

.914495 

Y4  Y5  . 

18418.419 

NOTE:  NSFS  SETS  OF  DATA  ARE  READ  HERE. 

NOTE:  MORE  THAN  ONE  CARD  MAY  BE  REQUIRED  FOR  XI  OR  YI . 

NOTE :  NXAPRX  VALUES  OF  X  AND  NF  VALUES  OF  Y  ARE  READ  FOR  EACH  SET  OF  DATA 

NOTE:  MANY  SIGNIFICANT  DIGITS  ARE  DESIRABLE. 


FIELD 

CONTENTS 

1-8 

XI: 

Ith  value  of  X,  I  =  1, NXAPRX. 

1-8 

YI: 

Ith  value  of  Y,  I  =  1,NF. 
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DATA  BLOCK  P  OMIT  IF  NSV  =  0  IN  BLOCK  B 
DESCRIPTION:  Sensitivity  objectives  (function  values) . 

FORMAT  AND  EXAMPLE 


1 

NSOBJ 

5 

2 

IPSENS 

0 

3 

4 

5  6  7 

8 

FORMAT 

2110 

NSNI 

3 

NSN2 

4 

NSN3 

5 

NSN4 

6 

NSN5 

7 

8110 

NOTE:  TWO  OR  MORE  CARDS  ARE  READ  HERE. 


FIELD  CONTENTS 

1  NSOBJ:  Number  of  separate  objective  functions  to  be  calculated  as 

functions  of  the  sensitivity  variables. 

2  IPSENS:  Print  control.  If  IPSENS.GT.O,  detailed  print  will  be  called 

at  each  step  in  the  sensitivity  analysis.  DEFAULT  =  No  print. 

1-8  NSNI:  Global  variable  number  associated  with  the  sensitivity 

objective  functions. 

REMARKS 

1)  More  than  eight  sensitivity  objectives  are  allowed.  Add  data  cards  as 
required  to  contain  data. 
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DATA  BLOCK 


Q  OMIT  IF  NSV  =  0  IN  BLOCK  B 


DESCRIPTION :  Sensitivity  variables . 

FORMAT  AND  EXAMPLE 


1 

I  SENS 

9 

2 

NSENS 

4 

3 

4  5 

6 

7 

8  FORMAT 

2110 

SNS1 

200. 

SNS2 

100. 

SNS3 

150. 

SNS4  .... 

250. 

. . . . 

.  8F10 

NOTE: 

READ  ONE 

SET  OF  DATA  FOR  EACH  OF  THE  NSV  SENSITIVITY  VARIABLES. 

NOTE: 

TWO  OR  MORE  CARDS  ARE  READ  FOR  EACH  SET  OF  DATA. 

FIELD 

CONTENTS 

1 

ISENS : 

Global  variable  number  associated  with  the  sensitivity 
variable . 

2 

NSENS: 

Number  of  values  of  this  sensitivity  variable  to  be  read  on 
the  next  card. 

1-8 

SENS I : 

Values  of  the  sensitivity  variable.  I  =  1, NSENS.  1=1 
correspond  to  the  nominal  value. 

REMARKS 

1)  More  than  eight  values  of  the  sensitivity  variable  are  allowed.  Add  data 
cards  as  required  to  contain  the  data. 


DATA  BLOCK  R  OMIT  IT  N2VAR  =  0  IN  BLOCK  B 
DESCRIPTION:  Two  variable  function  space  control  parameters. 

FORMAT  AND  EXAMPLE 

1  2  3  4  5  FORMAT 

N2VX  M2VX  N2VY  M2VY  IP2VAR  5110 


CM 

5  0 

FIELD 

CONTENTS 

1 

N2VX: 

Global  location 
function  space . 

of  the  X-variable  in 

the 

two  variable 

2 

M2VX: 

Number  of  values 

of  X  to  be  considered. 

3 

N2VY : 

Global  location 
function  space. 

of  the  Y-variable  in 

the 

two  variable 

4 

M2VY : 

Number  of  values 

of  Y  to  be  considered. 

5 


IP2VAR: 


Print  control.  If  IP2VAR.GT.0,  detailed  print  will  be  called 
at  each  step  (each  X-Y  combination).  DEFAULT  =  No  print. 


DATA  BLOCK 


S  OMIT  IF  N2VAR  =  0  IN  BLOCK  B 


DESCRIPTION:  Objective  functions  of  the  two  variable  function  space  study. 


FORMAT  AND  EXAMPLE 


1 

NZ1 

3 


2  3 

NZ2  NZ3 

4  5 


4  5  6  7  8  FORMAT 

NZ4  NZ5  .  8110 

6  7 


FIELD 


CONTENTS 


1-8  NZI : 


Global  location  corresponding  to  the  Ith  function  of  X  and  Y 
to  be  calculated.  N2VAR  values  are  read  here. 


REMARKS 

1)  More  than  eight  objective  functions  are  allowed.  Add  data  cards  as  required 
to  contain  the  data. 


DATA  BLOCK 


OMIT  IF  N2VAR  =  0  IN  BLOCK  B 


DESCRIPTION: 

Values  of  the  X-variable  in  a  two  variable  function  space  study. 

FORMAT  AND  EXAMPLE 

1 

2 

3 

45678  FORMAT 

XI 

X2 

X3 

X4  .  8F10 

0.5 

1.0 

1.5 

2.0 

FIELD 

CONTENTS 

1-8 

XI: 

Values  of  the  X-variable  in  the  two  variable  function  space 
M2VX  values  are  read  here. 

REMARKS 

1)  More  than  eight  values  are  allowed.  Add  data  cards  as  required  to  contain 
the  data. 
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DATA  BLOCK 


U  OMIT  IF  N2VAR  =  0  IN  BLOCK  B 


DESCRIPTION:  Values  of  the  Y-variable  in  a  two  variable  function  space  study. 

FORMAT  AND  EXAMPLE 


1 

2 

3  4  5  6  7 

8 

FORMAT 

Y1 

Y2 

Y3  Y4  Y5  . 

0  , 

8F10 

4.0 

8.0 

12.0  16.0  20.0 

FIELD 

CONTENTS 

1-8 

YI: 

Values  of  the  Y-variable  in  the 
N2VY  values  are  read  here. 

two 

variable  function  space 

REMARKS 

1)  More  than  eight  values  are  allowed.  Add  data  cards  as  required  to  contain 
the  data. 
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DATA  BLOCK 


V 


DESCRIPTION:  COPES  data  'END'  card. 

FORMAT  AND  EXAMPLE 

1  FORMAT 

END  3A1 

END 

FIELD  CONTENTS 

1  The  word  'END'  in  colums  1-3 

REMARKS 

1)  This  card  MUST  appear  at  the  end  of  the  COPES  data. 

2)  This  ends  the  COPES  input  data. 

3)  Data  for  the  user-supplied  routine,  ANALIZ,  follows  this 
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2.  TRANSONIC  ANALYSIS  INPUT  DATA 


NOTE: 

Excluding  literal  cards, 

all  input  data  cards  are  7F10.  format. 

CARD 

NUMBER 

CARD 

COLUMN 

VARIABLE 

NAME 

DESCRIPTION 

Card  1-A 

1-80 

TTLE 

Configuration  or  run  title  to  identify 
graphic  and  printed  output. 

Card  2-A 

1-10 

CASE 

=  1.  Isolated  Body,  used  for  input 

check  of  complex  body  definition. 
No  flow  solution.  (Omit  cards 
3-A,  4-A  and  all  cards  -C,  -W) . 


=  2.  Isolated  Wing  (omit  all  cards  -C, 
-B). 

=  3.  Wing-Body  (omit  cards  -C). 

=  4.  Isolated  wing-canard  (omit  all 
cards  -B) . 


=  5.  Wing-Body-Canard. 


11-20 

AMACH 

Mach  Number  (AMACH  <  1.0). 

21-30 

AOA 

Angle-of-Attack  (degrees). 

31-40 

WPO 

<-3.  Same  as  WPO  =  -2.  plus  omit  grid 

and  body  information  output. 

<-2.  Same  as  WPO  =  -1.  plus  omit  body 
0^  output. 

<-l.  Same  as  WPO  =  0.  plus  omit  Mach 
chart  output. 

<0.  No  crude  grid  output. 

<1.  Crude  grid  output  for  diagnostic 
purposes . 

>2.  Same  as  WPO  =  1.  plus  boundary 
layer  information. 

>3.  .  Same  as  WPO  =  2.  plus  print  C 
values  off  wing  and  canard.  ^ 
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CARD 

NUMBER 

Card  2-A 
(cont'd) 


Card  3-A 


NOTE: 


NOTE: 


* 


CARD  VARIABLE 

COLUMN  NAME  DESCRIPTION 


>4.  Same  as  WPO  =  3.  plus  print  solution 
convergence  information  for  every  span- 
wise  plane. 


41-50 

AXIT 

Number  of  initial  crude  grid  iterations. 

51-60 

AXITF 

Number  of  crude/fine  grid  iteration 
cycles . 

61-70 

VISMOD 

=1.  No  viscous  effects. 

=  2.  Viscous  effects  computed  at  end  of 
inviscid  analysis. 

=  3.  Inviscid/viscous  interaction. 

1-10 

FSAVE 

=  1.  Save  the  flow  solution  on  Unit  99. 

f  1.  Do  not  save  the  flow  solution. 

11-20 

FSTRT 

=  1.  Restart  the  flow  solution  from 

Unit  98. 

^  1.  Do  not  restart  the  flow  solution. 

See  Usage  Note  1. 

21-30 

CNVTST 

Convergence  test  based  on  average  flow 
solution  correction  (DPMAV).  Default 
is  1.0E-05. 

31-40 

FCASM2 

=  1.  Construct  wing-canard  grid,  but  run 
canard-off  solution.  Will  read 
wing-canard  of  wing-body-canard 
input  (CASE  =4.  or  5.  on  Card  2-A), 
construct  grids  and  then  do  wing 
alone  or  wing-body  solution  (CASE 
changed  to  CASE-2,  for  the  solution) 

$  1 .  Do  not  "turn  off"  canard. 

41-50 

PCTLE 

Crude  grid  leading  edge  spacing 
tolerance.  Default  is  0.01,  minimum 
allowed  is  0.005. 

See  Usage  Note  2. 
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1 


CARD 

NUMBER 

CARD 

COLUMN 

VARIABLE 

NAME 

DESCRIPTION 

Card  3-A 
(coat'd) 

51-60 

FCGRD 

=  1.  Places  the  canard  surface  in  the 

vertical  mid-location  of  the  canard 
fine  grid  system.  Appropriate  for 
lightly  or  negatively  loaded  canards. 

^1.  Places  the  canard  surface  at  the 
lower  quarter  vertical  location  in 
the  canard  fine  grid  system. 

61-70 

FMESH 

=  1.  Changes  crude  grid  dimensions  to 

51,21,27  X,Y,Z  points,  respectively. 

t  1.  Default  crude  grid  dimensions  of 

51,26,31  X,Y,Z  points,  respectively. 

Card  4-A 

1-10 

SREF 

Reference  area,  if  <0.,  code  will 
calculate. 

11-20 

AMAC 

Mean  aerodynamic  chord,  if  <0.,  code 
will  calculate. 

21-30 

ALAM 

Reference  taper  ratio,  if  <0.,  code  will 
calculate. 

31-40 

XMOM 

X-position  for  pitching  moment  reference. 

41-50 

ZMOM 

Z-position  for  pitching  moment  reference. 

51-60 

RE 

Reynolds  Number  xlO  based  on  AMAC. 

61-70 

FYINT 

=  0.  Nondimensional  ordinate  spanwise 
interpolation. 

=  1.  Physical  ordinate  spanwise 
interpolation. 

NOTE: 

See  Usage 

Note  3. 

NOTE: 

Omit  Card 

1-C  for  CASE<4. 

Card  1-C 

1-10 

ASECT 

Number  of  streamwise  sections  defining 
canard  planform  (2.<ASECT<20. ) . 

ASECT  for  canard  plus  ASECT  for  wing 
must  not  exceed  20. 

11-20 

AMIN 

Number  of  ordinates  defining  each 
canard  section  (ANIN<60). 
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CARD 

CARD 

VARIABLE 

NUMBER 

COLUMN 

NAME 

DESCRIPTION 

Card  1-C 
(cont'd) 

21-30 

ANOSW 

=  0.  Sharp  nose  canard  section. 

=  1.  Blunt  nose  canard  section. 

31-40 

ZWINGC 

Z-position  of  canard  (waterline). 

41-50 

XTRNC 

Transition  location  for  canard, 
streamwise. 

>  0.  fraction  of  chord 

<  0.  physical  distance  from  leading  edge 

=  0.  default  to  fixed  chord  fraction  of 
0.05. 

NOTE: 

See  Usage 

Note  4. 

51-60 

CINSDS 

Canard  incidence,  degrees. 

61-70 

SHIFTC 

Grid  shift  for  CASE>3. 

NOTE: 

See  Usage 

Note  5 . 

Card  1-W 

1-10 

ASECT 

Number  of  streamwise  sections  planform 
(2.<ASECT<20. ) . 

ASECT  for  canard  plus  ASECT  for  wing 
must  not  exceed  20. 

11-20 

ANIN 

Number  of  ordinates  defining  each  wing 
section  (ANIN  <60. ). 

21-30 

ANOSW 

=  0.  Sharp  nose  wing  sections. 

=  1.  Blunt  nose  wing  sections. 

31-40 

ZWING 

Z-position  of  wing  (waterline) . 

41-50 

XTRNV 

Transition  location  for  wing,  streamwise 

Same  usage  as  for  XTRNC  (See  Card  1-C 
and  Usage  Note  4.). 

51-60 

WINSDS 

Wing  incidence,  degrees. 

61-70 

SHIFTW 

Grid  shift  for  CASE<4. 

NOTE: 

See  Usage 

Note  5. 
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NOTE: 

Omit  Card 

Set  2-C  through  5 

-C  for  CASE<4.  Card  Set  2-C  through 

5-C  is  repeated  ASECT  (for 

canard)  times. 

Card  2-C 

1-10 

XPL 

Canard  section  leading  edge  (X-value) . 

11-20 

YP 

Canard  section  span  position  (Y-value). 
First  Y-value  must  be  0.0  (symmetry 
plane),  even  for  wing-body  case. 

21-30 

XPT 

Wing  section  trailing  edge  (X-value) . 

31-40 

TWIST 

Wing  section  local  incidence  (twist 
angle  in  degrees). 

41-50 

AKODE 

=  0.  Section  ordinates  identical  to 
preceding  section  (omit  cards  3-C 
through  5-C). 

=  1.  New  section  definition  expected  on 
cards  4-C  and  5-C. 

Card  3-C 

1-70 

XINWC 

Canard  section  x/c  coordinates  (cards 
3-C  defined  only  for  first  canard 
section,  ANIN  values). 

Card  4-C 

1-70 

YINU 

Canard  section  upper  surface  y/c 
coordinates  (ANIN  values). 

Card  5-C 

1-70 

YINL 

Canard  section  lower  surface  y/c 
coordinates  (ANIN  values). 

NOTE: 

Card  Set 

2-W  through  5-W  is 

repeated  ASECT  (for  wing)  times. 

Card  2-W 

1-10 

XPL 

Wing  section  leading  edge  (X-value). 

11-20 

YP 

Wing  section  span  position  (Y-value). 
First  Y-value  must  be  0.0  (symmetry 
plane),  even  for  wing-body  case. 

21-30 

XPT 

Wing  section  trailing  edge  (X-value) . 

31-40 

TWIST 

Wing  section  local  incidence  (twist 
angle  in  degrees). 

41-50 

AKODE 

=  Q.  Section  ordinates  identical  to 
preceding  section  (omit  cards  3-W 
through  5-W) 

=  1.  New  section  definition  expected  on 
cards  4-W  and  5-W. 
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CARD 

NUMBER 

CARD 

COLUMN 

VARIABLE 

NAME 

DESCRIPTION 

Card  3-W 

1-70 

XINW 

Wing  section  x/c  coordinates  (cards  3-W 
defined  only  for  first  wing  section, 
ANIN  values  expected). 

Card  4-W 

1-70 

YINU 

Wing  section  upper  surface  y/c  coordi¬ 
nates  (ANIN  values). 

Card  5-W 

1-70 

YINL 

Wing  section  lower  surface  y/c  coordi¬ 
nates  (ANIN  values). 

NOTE: 

Omit  card 

set  1-B  through 

13-B  for  CASE  =  2  or  CASE  =  4. 

Card  1-B 

1-10 

BKOD 

=  1.  Infinite  cylinder  (only  RADIUS 
need  be  input). 

=-l.  Same  as  BKOD  =1.  No  embedded  body 
grid.  Crude  grid  body  representa¬ 
tion  only. 


=  2.  Simple  axisymmetric  body  definition 
requested  (input  XB1N,  RIN  on  card(s) 
2-B  and  3-B. 

=-2.  Same  as  BKOD  =2.  No  embedded  body 
grid.  Crude  grid  body  representa¬ 
tion  only. 

=  3.  Complex  body  definition  requested 
(input  Quick-Geometry  model  on 
card(s)  4-B  through  13-B). 

=-3.  Same  as  BKOD  =3.  No  embedded  body 
grid.  Crude  grid  body  represen¬ 
tation  only. 


NOTE: 

BKOD>0 .  is 

available 

for  body  input  checkout  only.  Flow  solution  is 

available 

for  BK0D<0 . 

only. 

NOTE: 

See  Usage  Note  6. 

11-20 

BN0SE 

Body  nose  (X-value) 

for  BK0D=±2 .  or  ±3. 

21-30 

BTAIL 

Body  tail  (X-value) 

31-40 

BNIN 

Number  of  axisymmetric  body  coordinates 

to  be  input.  BNIN  <  60.  (for  BKOD 
=  ±2.  only). 

41-50 

RADIUS 

Cylinder  radius  for  BKOD  =  ±1.  only. 
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CARD 

CARD 

VARIABLE 

NUMBER 

COLUMN 

NAME 

DESCRIPTION 

Card  1-B 
(cont'd) 

51-60 

ANOSB 

=  0.  Sharp  nose  body 

BKOD  =  ±2.  only. 

=  1 .  Blunt  nose  body 

NOTE: 

Omit  card 

sets  2-B  and  3-B 

for  BKOD  =  ±1.  or  BKOD  =  ±3. 

Card(s)  2- 

-B  1-70 

XINB 

Axisymmetric  body  X-coordinates  (BNIN 
values) . 

Card(s)  3- 

-B  1-70 

RIN 

Axisymmetric  body  radii  (BNIN  values). 

NOTE: 

Omit  card 

sets  4-B  through 

13-B  for  BKOD  =  ±1.  or  BKOD  =  ±2. 

Card  4-B 

1-70 

VTITLE 

Quick-Geometry  model  title 

Card  5-B 

1-10 

ACSM 

Number  of  distinct  cross-section 
models  (ACSM  card  sets  6-B  and  7-B 
will  follow). 

Card  6-B 

1-10 

ADUM 

Running  count  of  current  cross-section 
model  (1-ACSM) . 

11-20 

AARC 

Number  of  arcs  in  current  cross-section 
model  (AARC  Card(s)  7-B  will  follow). 

21-60 

CTITLE 

Title  or  descriptor  of  current  cross- 
section  model. 

Card  7-B 

1-8 

ARCNAM 

Arc  or  component  name. 

11-14 

ASHAPE 

Arc  or  component  shape. 

21-28 

PNTNAM(l) 

Control  point  name  for  beginning  of 
this  arc. 

32-28 

PNTNAM(2) 

Control  point  name  for  termination  of 
this  arc. 

41-48 

PNTNAM(3) 

Slope  control  point  name  for  this  arc, 
if  required. 

Card  8-B 

1-10 

ANTCSM 

Number  of  cross-section  models  to  define 
entire  body  (ANTCSM  card(s)  9-B  will 
follow) . 

1 


CARD 

NUMBER 


CARD 

COLUMN 


VARIABLE 

NAME 


DESCRIPTION 


Card 

9-B 

1-10 

ADUM 

Running  count  of  current  cross-section 
model  (1-ANTCSM) . 

11-20 

AMODEL 

Index  corresponding  to  already  defined 
cross-section  models  (between  1  and 
ACSM) . 

21-30 

XCSMS1 

Starting  X-station  for  current  cross- 
section  model. 

31-40 

XCSMS2 

Ending  X-station  for  current  cross- 
section  model . 

Card 

10-B 

1-10 

BLINE 

Number  of  body  line  models  to  be  defined 
by  segments  (BLINE  card  set  pairs  11-B 
and  12-B  follow) . 

11-20 

ALIAS 

Number  of  body  line  models  to  be  aliased 
(Input  ALIAS  card(s)  13-B  below). 

Card 

11-B 

1-10 

BLSEG 

Number  of  segment(s)  defining  body 
line  model . 

11 

BYORZ 

The  letter  Y  or  Z  indicates  which  data 
definition  is  to  follow. 

12-19 

BNAME 

Body  line  name  to  be  defined. 

Card 

12-B 

1-4 

SSHAPE 

Segment  shape. 

11-20 

D(l) 

X-station  for  beginning  of  segment. 

21-30 

D(2) 

Y  or  Z  value  corresponding  to  D(l). 

31-40 

D  ( 3 ) 

X-station  for  termination  of  segment. 

41-50 

D(4) 

Y  or  Z  value  corresponding  to  D(3). 

51-60 

D(5) 

X-station  for  segment  slope  control 
point. 

61-70 

D(6) 

Y  or  Z  value  corresponding  to  D(5). 

CARD 

NUMBER 

CARD 

COLUMN 

VARIABLE 

NAME 

DESCRIPTION 

Card  13-B 

11 

BYORZ 

The  letter  Y  or  Z  indicates  which  data 
definition  is  to  follow. 

12-19 

BNAME 

Body  line  name  to  be  defined. 

21 

AYORZ 

The  letter  Y  or  Z  indicates  which 
definition  is  to  be  used  for  aliasing. 

22-29 

ANAME 

Body  line  name  to  which  BNAME  is 
aliased. 

USAGE  NOTES 


1 .  Saved  Solution 

When  an  analysis  consists  of  both  crude  only  and  then  crude/fine 
iterations,  a  saved  solution  is  written  on  unit  99  immediately  before  the 
crude/fine  iterations  begin.  At  the  end  of  the  crude/fine  iterations,  the 
previous  saved  solution  is  overwritten  with  the  most  recent  results.  In  this 
way,  an  abnormal  termination  of  the  crude/fine  iterations  will  have  the  crude 
iteration  results  saved  on  unit  99.  With  an  abnormal  termination, 
appropriate  job  control  cards  may  be  needed  to  "permanently"  save  the  unit  99 
data . 

2.  PCTLE  Parameter 


The  code  will  shift  the  input  geometry  to  find  a  streamwise  location  in 
the  crude  grid  where  the  first  points  at  each  span  station  on  the  wing  and 
canard  is  not  "too  close"  to  the  leading  edge.  (Riegel's  rule  is  not  used  for 
the  boundary  conditions.)  The  PCTLE  parameter  is  the  required  minimum 
distance  from  the  leading  edge.  See  also  Usage  Note  5. 
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3. 


FYINT  Parameter 


The  input  wing  and  canard  section  ordinates  (YINU  AND  YINL)  are  lin¬ 
early  interpolated  from  the  input  span  stations  to  the  analysis  span  sta¬ 
tions.  The  method  of  interpolation  is  controlled  by  the  parameter  FYINT. 
If  FYINT  =  0.  nondimens ional  spanwise  interpolation  of  the  ordinates  is 
used: 


(z/c)Y  =  (z/c)Y1(l .-R)  +  (z/c)Y2(R) 

If  FYINT  =  1.  physical  spanwise  interpolation  of  the  ordinates  is  used: 
(z/c)y  =  (z/c)Y1(l ,-R) (cY1/cY)  +  (z/c)y2(R)(cY2/cY) 

where 

c  =  local  chord  Y  =  interpolated  span  station 

R  =  (Y-Yl)/ (Y2-Y1)  Y1,Y2  =  input  span  stations 

z/c  =  nondimensional  section  ordinate 

The  first  formula  is  the  more  usual  analysis  code  interpolation,  while 
the  second  formula  corresponds  to  manufacturing  lofting  methods.  Waggoner 
[10]  has  shown  the  difference  to  be  significant  for  his  application.  The 
difference  between  the  two  formulas  becomes  greater  for  more  highly  tapered 
wings  and  fewer  input  span  stations. 

4.  Boundary  Layer  Transition 

The  fraction  of  chord  designation  results  in  transition  specified  at 
the  same  percentage  chord  location  at  all  span  stations.  Thus,  the  physical 
distance  from  the  leading  edge  decreases  as  the  local  chord  length  de¬ 
creases.  The  physical  distance  designation  results  in  transition  specified 
at  the  same  fixed  distance  from  the  leading  edge.  Thus,  the  percentage 
chord  from  the  leading  edge  increases  as  the  local  chord  length  decreases. 

5.  Grid  Shift  Parameters 


Variables  SHIFTC  and  SHIFTW  can  be  used  to  set  the  initial  streamwise 
placement  of  the  configuration.  (The  grid  is  fixed  in  location  and  the 
input  geometry  is  shifted.)  The  code  will  determinate  the  appropriate 
shift  with  a  search  procedure.  After  the  first  analysis  run,  the  resulting 
shift  parameter  can  be  input  to  eliminate  the  search  in  subsequent  runs. 

The  shift  parameter  can  also  be  used  to  evaluate  different  grid  placements. 
If  the  grid  placement  search  should  fail,  the  shift  parameter  can  be  used 
to  start  the  search  at  different  locations,  possibly  finding  a  satisfactory 
grid  placement. 
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uum **  * 


Body/Fuselage  Geometry  Model 

The  present  method  allows  complex  three-dimensional  geometries  to  be 
input,  processed  and  converted  into  a  suitable  array  of  boundary  conditions 
for  analysis.  Although  the  input  or  modeling  of  complex  body  shapes  is 
extremely  error  prone  and  certain  applications  might  not  warrant  this  level 
of  effort,  it  is  necessary  in  aircraft  applications  when  fuselage  contours 
(e.g.,  canopies,  fairings)  are  required.  The  code  does  not  provide  an 
embedded  fine  grid  body  analysis.  The  embedded  body  grid  option  is  pro¬ 
vided  to  allow  a  detailed  body  input  definition  checkout.  This  will  be 
most  useful  if  the  graphical  output  is  available.  The  following  discussion 
is  excerpted  from  Reference  [ 1 ] .  Additional  information  and  examples  will 
be  found  in  the  reference. 

The  complex  fuselage  modeling  system  has  been  named  "Quick- Geometry" 
by  its  developers,  Vachris  and  Yaeger  [6].  A  detailed  User's  Guide  for  the 
Quick-Geometry  system  can  be  found  in  the  appendix  of  Reference  [7].  This 
system  was  originally  developed  for  the  geometric  modeling  of  wing-body 
shapes.  Since  only  fuselage  shapes  are  of  concern  here,  many  of  the  more 
sophisticated  options  including  fillets  and  patches  will  not  be  described 
in  the  paragraphs  which  follow.  In  addition,  if  References  [6]  and  [7]  are 
being  used  to  augment  the  modeling  description  provided  herein,  it  should 
be  noted  that  the  input  format  has  been  modified  to  be  more  consistent  with 
that  of  the  basic  transonic  wing-body  code. 

The  geometry  package  requires  that  certain  body  lines  and  cross- 
section  lines  be  defined.  The  body  lines  and  cross-section  lines  may  be 
likened  to  the  stringers  and  bulkheads,  respectively,  used  in  fuselage 
construction.  These  line  models  are  defined  by  a  combination  of  simple 
curves  (i.e.,  lines,  ellipses,  cubics).  They  are  taken  together  to  provide 
a  continuous  analytical  model  of  the  surface  geometry.  Slopes  and  normals 
are  developed  analytically.  Either  discontinuous  intersections  or  smooth 
fairings  can  be  modeled  and  enforced. 

Two  different  coordinate  systems  are  employed.  Geometry  definition  is 
performed  in  a  Cartesian  coordinate  system  (x,  y,  z) ,  while  interrogation 
of  the  model  for  body  boundary  conditions  is  performed  in  cylindrical 
coordinates  (X,  R,  0).  This  results  in  the  use  of  a  plane  of  symmetry  map 
axis,  the  height  of  which  usually  corresponds  to  the  position  of  the  max- 
half-breadth  line.  It  is  required  that  the  configuration  radius  at  any 
cross-sectional  cut  be  a  single  valued  function  of  the  angle  0.  These 
definition  lines  and  coordinate  systems  are  illustrated  in  Figure  6. 

A  minimum  of  four  body  lines  are  required  for  the  simplest  fuselage. 
These  are:  1)  top  centerline,  2)  bottom  centerline,  3)  max-half-breadth 
line,  and  4)  the  map  axis.  Each  body  line  must  be  defined  by  both  its  Y 
and  Z  values  over  the  full  range  of  X  (between  fuselage  nose  and  tail). 
Similarly,  a  minimum  of  two  cross-section  line  segments  are  required  for 
each  different  cross-section  line  model.  These  are  1)  body  upper,  and  2) 
body  lower. 
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Both  body  lines  and  cross-section  lines  are  specified  by  defining  key 
arc  or  segment  shapes  and  their  accompanying  limiters.  The  segment  shape 
boundary  conditions  that  are  used  to  determine  the  coefficients  of  the 
slope  equation  are  the  1)  origin  point,  2)  termination  point,  and  3)  slope 
control  point.  The  slope  control  point  lies  at  the  intersection  of  the 
line  which  is  tangent  to  the  segment  shape  at  the  origin  point  and  the  line 
which  is  tangent  to  the  segment  shape  at  the  termination  point  (see  Figure 
7).  The  slope  control  point  is  a  very  convenient  way  of  specifying  slope 
conditions.  In  particular,  it  allows  for  the  simultaneous  specification  of 
slope  conditions  at  both  ends  of  the  segment. 

Figure  7  is  a  schematic  illustrating  the  component  build-up  of  a 
particular  body  line  and  cross-section  line  model.  Naturally,  LINE  seg¬ 
ments  do  not  require  a  slope  control  point.  In  this  case,  the  portion  of 
the  body  top  center  line  illustrated  requires  four  body  line  segments  and 
the  cross-section  is  constructed  with  two  arcs  (two  is  the  minimum  number 
allowed) . 

The  arc  shapes  used  for  defining  a  cross-section  line  model  are  listed 
in  Table  1. 


TABLE  1.  CROSS-SECTION  ARC  SHAPES 


SHAPE 

KEYWORD 

EQUATION 

Line 

LINE 

Ay 

+  Bz  +  C  =  0 

Ellipse 

(Concave  to  origin) 

ELLI 

(y 

•  V2  .  (Z  ‘  Z0} 

2  +  2 

AZ  BZ 

Ellipse 

(Convex  to  origin) 

ELLO 

Same  as  ELLI . 

They  are  input  in  an  order  which  starts  at  the  body  bottom  centerline  and 
proceeds  to  the  body  top  center  line.  The  segment  shapes  used  for  defining  a 
body  line  model  are  listed  in  Table  2. 
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TABLE  2.  BODY  LINE  SEGMENTS 


SHAPE  KEYWORD 


EQUATION 


Line 

LINE 

X-PARABOLA 

XPAR 

Y-PARABOLA 

YPAR 

X-ELLIPSE 

ELLX 

Y-ELLIPSE 

ELLY 

Ax  +  By  =  0 

Ax  +  By  +  y2  =  0 

Ax  +  By  +  x2  =  0 

Ax  +  By  +  Cx2  +  y2  =  0 

Ax  +  By  +  Cy2  +  x2  =  0 


CUBIC  CUBI  Ax  +  By  +  Cx2  +  x3  =  0 


Cross-sections  arcs  are  input  in  their  order  of  appearance.  However,  body 
line  segments  are  defined  along  with  an  index  which  establishes  their  order  in  the 
X-direction.  In  addition,  body  lines  may  be  aliased  to  other  body  lines.  This 
allows  two  body  line  models  to  have  identical  mathematical  representations  without 
repeating  the  body  line  segment  input.  For  example,  the  Z-value  of  the  map  axis 
(ZMAP)  is  typically  the  Z-value  of  the  maximum  half-breadth  (ZMHB) .  The  two  are 
made  identical  by  aliasing  the  two  body  line  model  names  ZMAP  and  ZMHB. 


It  should  be  noted  that  cross-sections  are  defined  only  in  terms  of  named 
component  arcs  (arc  shape  table)  and  named  control  points.  On  the  other  hand, 
body  lines  are  defined  mathematically  by  coordinates  over  the  length  of  the 
configuration  for  which  they  are  required.  At  a  given  X-station,  the  body  lines 
are  interrogated  to  give  the  key  control  points  required  to  construct  the  cross- 
sectional  arcs. 
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SECTION  V 


CODE  OPERATION 

The  program  reads  data  from  unit  5  and  writes  output  on  unit  6.  Units  20  and 
40  are  used  as  scratch  files  in  the  COPES  routines.  The  scratch  file  numbers  may 
be  changed  by  changing  two  cards  at  the  beginning  of  the  COPES  program.  The 
analysis  code  stores  data  on  unit  numbers  1,7,8,15,85,98  and  99.  Unit  15  is  used 
by  the  Ames  plotting  software  for  graphical  output.  It  has  been  left  available  to 
facilitate  incorporation  of  graphic  capability  by  other  users. 

The  computer  code  is  written  in  FORTRAN,  employing  the  CDC  overlay  structure. 
Storage  requirements  on  a  CDC  7600  are  151000-  small  core  memory  and  770000-  large 
core  memory.  Typical  CPU  times  on  the  7600  for  a  wing-body-canard  analysis  with 
viscous  effects  is  eleven  minutes.  The  CPU  time  for  an  optimization  run  will  vary 
according  to  the  complexity  of  the  optimization  problem,  as  discussed  in 
Section  3. 

The  overlay  structure  is  described  in  Table  3.  The  overlay  sizes  are  shown 
in  Table  4.  The  job  control  language  (JCL)  for  operation  on  the  Ames  CDC  7600  is 
shown  in  Figure  8.  The  job  card  and  account  card  will  vary  with  the  individual 
user.  Compilation  of  the  entire  program  requires  one  minutes  of  CPU  time.  To  save 
CPU  time  for  short  test  runs,  a  binary  file  of  the  code  is  stored  as  a  permanent 
file.  Thus,  the  JCL  uses  an  UPDATE, Q  command  followed  by  FTN  and  C0PYL  to  change, 
compile  and  replace  only  those  routines  being  modified  for  a  particular  run.  The 
numerical  optimization  will  require  changes  to  the  subroutines  OGEOM,  MOD  and 
ANALIZ  (all  in  Overlay  (0,0)).  For  simple  analysis,  no  changes  need  be  made, 
except  perhaps  to  BLOCK  DATA  to  change  some  defaults.  The  UPDATE, Q  always 
requires  some  input.  The  *C0MPILE  INTEG  card  in  Figure  8  is  shown  as  an  example 
that  satisfies  this  requirement. 


TABLE  3.  OVERLAY  STRUCTURE 


OVERLAY  CONTENTS 

(0,0)  COPES  control  program,  subroutines  needed  throughout  the  code 


(1,0) 

(1,1) 

(1,2) 


Input  control  program,  QUICK  geometry  look-up  routines 
Input  processing 

Input  geometry  plotting,  body  area  calculations 


(2,0) 

(2,1) 

(2,2) 

(2,3) 


Transonic  analysis  control  program,  viscous  calculation  routines 

Potential  flow  solution  routines 

Pressure,  force  and  moment  calculation  and  output 

Output  plotting  program 


(3,0) 


Previously  used  for  opening  and  closing  plotting  files.  Currently 
not  used. 


(4,0)  Body  force  and  moment  calculation 

(5,0)  C0NMIN  optimization  routines 
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TABLE  4.  OVERLAY  SIZES 


OVERLAY 

SMALL  CORE  SIZE 

LARGE  CORE  SIZE 

(0,0) 

66506 

766710 

(1,0) 

101300 

766710 

(1,1) 

125321 

767010 

(1,2) 

104642 

766710 

(2,0) 

105477 

766710 

(2,1) 

137421 

766710 

(2,2) 

136514 

766710 

(2,3) 

112151 

766710 

(3,0) 

66517 

766710 

(4,0) 

101261 

766710 

(5,0) 

100524 

766710 

The  Quick-Geometry  routines  are  needed  for  the  input  geometry  processing  and 
the  body  force  and  moment  calculations.  The  small  core  restrictions  do  not  allow 
all  the  routines  10  be  in  0verlay(0,0) .  The  COPYBR  cards  are  used  to  copy  the 
necessary  routiner:  into  the  overlays  where  needed  (Overlays  (1,0),  (2,2)  and 
(4,0)). 

Not  shown  in  the  JC'L  are  the  cards  for  permanently  saving  the  flow  solution. 
The  Applied  Computational  Aerodynamics  Branch  at  Ames  can  provide  the  necessary 
authorization  and  information. 

If  the  code  is  converted  for  use  on  IBM  type  computers,  several  Quick- 
Geometry  variables  need  special  treatment.  These  variables  and  the  sub-routines 
where  they  appear  are  listed  in  Table  5.  The  REAL*8  declaration  is  needed  to  have 
sufficient  word  length  for  line  model  labeling  and  increase  the  accuracy  of  the 
body  model  calculations. 
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TABLE  5.  REAL*8  VARIABLES  FOR  IBM  USAGE  OF  QUICK- GEOMETRY 


SUBROUTINE  NAME 

VARIABLE  NAME 

CURVES 

A,  B,  C,  Y,  T,  X,  FACT,  RFACT,  S 

MODTV 

SUM,  ONE 

VDOTZ 

C,  ONE 

QWIKDE 

CPNTNM,  COMPNM 

BLMCHK 

CPNTNM,  COMPNM,  ANAME,  BNAME,  BLMNAM,  EQUNB,  EQUNA 

BLMDEF 

CPNTNM,  COMPNM,  ANAME,  BNAME,  BLMNAM,  BLANK2 

CSMCHK 

CPNTNM,  COMPNM,  ARCNAM,  PNTNAM , 
ARCNM 

BLANK2,  ZMAPNM, 

CSMDEF 

CPNTNM,  COMPNM,  ARCNAM,  PNTNAM, 
ARCNM,  AMAPAX 

BLANK2 ,  ZMAPNM, 

GEMOUT 

CPNTNM,  COMPNM,  ARCNM,  PNTNAM, 
ARCNM 

BLANK2,  AMAPNM, 

DSETUP 

ALABLE,  BLABLE 

The  program 
shown  in  Figure  9 

output  listing  of  the  input  data  for  a  sample 
.  In  front  of  this  input  were  the  COPES  control 

analysis  run  is 
cards : 

TITLE 

1. 

END 


with  no  delimiters  in  between.  The  input  example  uses  NACA0010  wing  sections. 
Following  the  input  data  listing  is  a  data  read  echo  showing  the  information  the 
program  has  read  from  the  input  cards  (Figure  10).  Portions  of  the  grid  gene¬ 
ration  output  are  shown  in  Figures  11  and  12. 

In  this  example,  the  flowfield  solution  begins  with  five  crude  grid  itera¬ 
tions,  and  then  continues  with  five  fine  grid/crude  grid  iterations.  The  program 
output  for  this  is  shown  in  Figure  13.  When  the  flow  solution  is  completed, 
detailed  flow  information  is  output  at  each  analysis  station  on  the  wing  and 
canard.  An  example  for  station  11  is  shown  in  Figure  14. 

Force  and  moment  output  for  the  wing,  canard  and  body  is  shown  in  Figure  15. 
Figure  16  shows  the  final  wing-body-canard  spanload  output.  All  of  the  results 
shown  in  Figure  10  through  16  were  for  the  analysis  input  data  of  Figure  9. 
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1 


AGn 

A1 

A2 

A3 

CCD/CAV 

CCF/CAV 

CCL/CAV 

CCM/CAV/MAC 

CD 

CDINT 

CF 

CIR 

CL 

CLINT 

CM 

CMLOC 

CP 

Cl 

C3 

DPM 

DPMAV 

DRAG 

D1 

D3 


APPENDIX  -  GLOSSARY  OF  OUTPUT  VARIABLES 

n=l,2,3,4,5  Coefficients  aQ  in  equation  (7). 

Coefficient  A^^  in  equation  (6). 

Coefficient  in  equation  (8) . 

Coefficient  A^  in  equation  (9). 

Spanload  drag. 

Spanload  skin  friction. 

Spanload  lift. 

Spanload  pitching  moment. 

Drag  coefficient  (local  or  total,  as  appropriate). 

Local  drag  coefficient  in  spanload  tables. 

Skin  friction  coefficient. 

Circulation. 

Lift  coefficient  (local  or  total  as  appropriate). 

Local  lift  coefficient  in  spanload  tables. 

Pitching  moment  coefficient. 

Local  pitching  moment  coefficient  (about  quarter  chord  point). 
Pressure  coefficient. 

Coefficient  C^  in  equation  (8). 

Coefficient  C^  in  equation  (8). 

Maximum  correction  to  flow  potential  for  the  current  iteration. 

Average  of  the  absolute  value  of  all  the  flowfield  corrections  for 
the  current  iteration. 

Local  drag  contribution  in  the  body  force  output. 

Coefficient  D^  in  equation  (9). 

Coefficient  D^  in  equation  (9). 


.  1 

!  \ 


i 

3 


1 


H 


HOCBUfiB  PAOK  BUNK-MOT  FILMED 


ETA 

I 

IL 

INOSE 

INOSEC 

IT 

ITAIL 

ITAILC 

ITER 

J 

JSD 

KLO 

LOAD 

NSP 

PCT 

PHI 

RMAX 

RSDAV 

SEXPC 

SHIFT 

U 

V 

WCORD 

X 

XI 


Spanwise  computational  ordinate. 

General  streamwise  index. 

Streamwise  index  of  leading  edge  point  in  crude  grid. 

Streamwise  index  of  the  first  point  on  the  body  in  the  body  fine 
grid. 

Streamwise  index  of  the  first  point  on  the  body  in  the  crude  grid. 

Streamwise  index  of  the  trailing  edge  point  at  an  analysis  station 
in  the  crude  grid. 

Streamwise  index  of  the  last  point  on  the  body  in  the  body  fine 
grid. 

Streamwise  index  of  the  last  point  on  the  body  in  the  crude  grid. 
Iteration  count. 

General  spanwise  index. 

Spanwise  index  of  first  solution  plane  outside  the  body 
computational  surface. 

Vertical  index  for  bottom  of  body  computational  surface. 

Local  lift  contribution  in  body  force  output. 

Number  of  supersonic  points  in  the  crude  grid. 

Fractional  cord  distance  from  leading  or  trailing  edge  during 
crude  grid  shift  search. 

Pertubation  velocity  potential. 

Maximum  body  radius . 

Average  of  the  absolute  value  of  the  error  in  the  finitie 
difference  equation  at  all  grid  points  for  the  current  iteration. 

Canard  exposed  area. 

Grid  shift  parameter. 

Streamwise  pertubation  velocity. 

Spanwise  pertubation  velocity. 

Local  cord  for  wing  or  canard. 

Streamwise  physical  coordinate. 

Streamwise  computational  coordinate. 
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XLE 

XNOSE 

XTA1L 

XTE 

XWF 

X/C 

Y 

Y/C 

Z 

ZETA 

2Y/B 


Wing  or  canard  leading  edge  X- location. 

Body  nose  X-location. 

Body  tail  X-location. 

Wing  or  canard  trailing  edge  location. 

Wing  or  canard  fine  grid  streamwise  physical  coordinate 
Percent  fraction  of  local  cord. 

Spanwise  physical  coordinate. 

Percent  chord  wing  or  canard  section  ordinate. 

Vertical  physical  coordinate. 

Vertical  computational  coordinate. 

Fraction  of  semispan  (on  wing  or  canard  as  appropriate) 
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Figure  2.  Multiple  wake  vertical  line  relaxation  regions,  crude  grid 
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Figure  5.  Fine/crude  grid  interface  (from  [1]) 
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Figure  6.  i)uick  geometry  model  lines  and  coordinate  system  (from  [1]  ) 
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Figure  7.  ljuick  geometry  body  and  cross-section  line  models  (from  [1]) 
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Figure  8  -  Job  control  language 
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Figure  13.  Analysis  code  data  read  echo 
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Figure  11.  Crude  grid  generation  output 
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Figure  11.  (Continued) 
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Figure  12.  Example  of  fine  grid  output 
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Figure  13.  (Concluded) 
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Figure  14.  Example  of  fine  grid  solution  output 
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Figure  14.  (Concluded) 


WING  EXPOSED  AREA  =  20361.104 

WING  REFERENCE  AREA  =  30237.000 

WING  ASPECT  RATIO  =3.56 

WING  TAPER  RATIO  =  .33 

W I NG  MEAN  AERODYN AM  I C  CHORD  =101.375 

WING  AVERAGE  CHORD  =  92.186 

X-POS I T I  ON  FOR  PITCHING  MOMENT  =  .000 

Z-POSITION  FOR  PITCHING  MOMENT  =  0.000 

WING  FORCE  AND  MOMENT  COEFFICIENTS 

WING  CL  WING  CM  WING  CD 
.04693  -.21111  .00061 


WING  PRESSURE  DRAG  =  .00061 

WING  FRICTION  DRAG  =  0  00000 

PITCHING  MOMENT  DUE  TO  DRAG  =  -.00018 

CANARD  FORCE  AND  MOMENT  BASED  ON  SEXPC  =  5431.315 

CL  CM  CD (PRESSURE)  CD (FRICTION) 

-.02768  .10143  -.00371  0.00000 

CANARD  FORCE  AND  MOMENT  COEFFICIENTS 

CANARD  CL  CANARD  CM  CANARD  CD 

-.00497  .01822  -.00067 


Figure  15.  Force  and  moment  output 
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CANARD  PRESSURE  DRAG  =  -.00067 

CANARD  FRICTION  DRAG  =  0.00000 

PITCHING  MOMENT  DUE  TO  DRAG  =  0.00000 

TOTAL  FORCE  AND  MOMENT  COEFFICIENTS 

TOTAL  CL  TOTAL  CM  TOTAL  CD 
.04196  -.19290  -.00006 


TOTAL  PRESSURE  DRAG  -  -.00006 
TOTAL  FRICTION  DRAG  =  0.00000 


Figure  15,  (Continued) 
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BODY  LENGTH  *  600.000 

BODY  WETTED  AREA  -  93000.300 
BODY  PROJECTED  AREA  =  29459.057 
BODY  MAX.  CROSS-SECTIONAL  AREA  =  4300.840 

X-POSITION  ABOUT  WHICH  MOMENTS  ARE  COMPUTED 


BODY  FORCE  AND  MOMENT  COEFFICIENTS 


BODY  CL 
.00819 


BODY  CM  BODY  CD 

-.03018  .00019 


BODY  PRESSURE  DRAG  «  .00019 

BODY  FRICTION  DRAG  »  0.00000 


Figure  15 


.  000 


.  (Concluded) 


4 
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TOTAL  CONFIGURATION  COEFFICIENTS 
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Figure  16.  Spanload  output 


